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We study the equilibrium statics and nonequilibrium driven dynamics of flux line liquids in pres- 
ence of a random pinning potential. Under the assumption of replica symmetry, we find in the static 
case using a replica Gaussian variational method that the only effect of disorder is to increase the tilt 
modulus and the confining "mass" of the internal modes of the flux lines, thus decreasing their ther- 
mal wandering. In the nonequilibrium, driven case, we derive the long scale, coarse-grained equation 
of motion of the vortices in presence of disorder, which apart from new Kardar-Parisi-Zhang non- 
linearities, has the same form as the equation of motion for unpinned vortices, with renormalized 
coefficients. This implies, in particular, that the structure factor of a disordered vortex liquid has 
the same functional form as in the absence of pinning, in disagreement with the results of previous 
hydrodynamic methods. The expression of the static structure factor derived within our approach is 
consistent both with experimental data and with the standard theory of elasticity of vortex lattices. 

PACS numbers: 74.20.De, 74.25.Qt 
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I. INTRODUCTION 

During the past fifteen years, the study of the proper- 
ties of flux line liquids in high temperature superconduc- 
tors (HTSCs) has been one of the most active areas of 
research in vortex phenomenology. Yet, despite an im- 
pressive body of literature, both experimental and the- 
oretical, and despite a relatively good understanding of 
the overall behavior and macroscopic properties of vor- 
tex liquids, it seems that the important question of the 
actual microscopic correlations of flux line trajectories 
inside such liquids has not been fully understood yet. 
Indeed, of the several theoretical approaches that have 
been used to study the properties of liquid vortex matter, 
one particular approach, which has had a rather strong 
impact on our present understanding of flux liquids in 
HTSCs, is the boson mapping, developed by several au- 
thors, 0, 0, S QL which is based on the observation 
that there is a formal mapping between the parti- 
tion function of a three-dimensional system of interact- 
ing flux lines, and the imaginary-time partition function 
of quantum bosons in two-dimensions. Although the bo- 
son mapping is ultimately used to find density-density 
correlation functions, and does not contain, after coarse- 
graining, any detailed information about flux-line trajec- 
tories, it has been argued, based on the behavior of 
the structure factor derived in this and other hydrody- 
namic approaches, 8] that flux lines wander throughout 
the sample in a random-walk-like fashion, much in the 
same way as in a hypothetical "ideal gas" of noninteract- 
ing vortices. This implies, in particular, that the internal 
fluctuations of flux lines have an average spatial extent 
which diverges with the sample thickness L. 

The above interpretation, and in fact the whole hydro- 
dynamic approach to flux line liquids, suffers from a num- 
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ber of inconsistencies which have been pointed out and 
discussed in detail in two recent papers by the author. 
0, In these two articles, the author has proposed a 
new approach to study three-dimensional flux-line liquids 
in type II superconductors which, instead of the density, 
uses the actual conformation variables of vortices as the 
fundamental dynamical variables of the flux-line system. 
This new approach, which makes contact with the stan- 
dard theory of classical fluids, is based on the separation 
of dynamical variables of flux lines into center of mass 
(cm. for short) and internal modes, and on the obser- 
vation that the repulsive interactions between flux lines 
must lead to a certain degree of confinement of the inter- 
nal modes, whose fluctuations are shown to be bounded 
and no- longer diverge with the sample thickness. 0, IT(tI | 
This picture is obviously in contradiction with the results 
of the boson mapping of refs. 0, and 0. 

In this paper we wish to generalize the method olog y 
developed in these previous studies, (refs. [tj and |10|). 
to study the statics and dynamics of vortex liquids in 
presence of a random pinning potential. In the static 
case, we shall show in particular that disorder leads to 
an enhancement of the tilt modulus and the confining 
"mass" of the internal modes of flux lines, thereby re- 
ducing their thermal wandering, but otherwise leaves the 
analytic form of the structure factor S(r, z) unchanged, 
contrarily to what has been argued by previous authors 
based on boson mapping methods 0,0 or other hydrody- 
namic approximations. 8] These static results are then 
generalized to the dynamics of driven vortex liquids in 
the presence of a random pinning potential. In contrast 
to earlier studies of this system, in which the density was 
used as the fundamental dynamical variable, in our ap- 
proach flux line trajectories, which are the true dynam- 
ical variables of the system, are used throughout. This 
enables us to derive the coarse-grained, large scale equa- 
tion of motion of vortices in the flux liquid in the presence 
of disorder, in close analogy with earlier work on driven, 
disordered flux lattices, fill 111 111 
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This article is organized as follows. In Sec. [HJ we use 
a simple Larkin analysis to study the effect of disorder 
on vortex liquids. This Larkin analysis is then refined in 
Sec. IIIII where we use the replica Gaussian variational 
method |l5| to properly include the effect of the rela- 
tively strong thermal fluctuations which are a common 
characteristic of HTSCs. In Sec. |lVl we construct an ac- 
tion formulation for the dynamics of flux lines in a vortex 
liquid. Then, before considering the nonequilibrium case 
of a driven disordered flux liquid, we shall first be inter- 
ested in the equilibrium dynamics of flux liquids, which 
we will investigate in quite some detail in Sec. V. Such 
an investigation is not only a natural step toward the 
more complex disordered case, but is also necessary for 
the developments to follow, since one needs to correctly 
specify the near-equilibrium dynamics of the unpinned 
interacting liquid in order to be able to tackle the out- 
of-equilibrium driven, disordered case. In Sec. IVII we 
derive the coarse-grained dynamics of driven, disordered 
flux liquids in the limit of high drives before deriving the 
structure factor of pinned flux liquids in Sec. IV 111 Sec. 
IVIIII contains a discussion of our results along with our 
conclusions. 



II. FLUX LINE LIQUID IN THE PRESENCE OF 
DISORDER: PERTURB ATIVE ANALYSIS 



We thus consider a flux-line liquid in d = d± + 1 di- 
mensions (we use d± = 2 in all explicit calculations), in 
presence of an external pinning potential. Our starting 
point is the Hamiltonian: P, |"~ 



"■l/^dME^w-iM)} 

1=1 

N 

+ E J dzV d (Ti(z),z) , (1) 



i=l ■ 

where the e?x-dimcnsional vector Vi{z) parametrizes the 
trajectory of the i-th flux-line as it traverses the super- 
conducting sample, K is the tilt modulus of the flux lines, 
V{r) = 2eoK (r / A) is the interaction potential between 
flux line elements at equal height, and Vd(r, z) is a ran- 
dom pinning potential. In the above expression of V(r), 
A is the London penetration de pth in the (ab) planes, Kq 
is a modified Bessel function, and £n = ((/>o/47rA) 2 , 
where 4>q = hc/2e is the flux quantum. |17| In equation 
(Q) and all equations below, the origin of heights is taken 
to be at the center of the sample, and all z integrals are 
taken from — L/2 to L/2 (L is the sample thickness). For 
simplicity, we shall consider that the probability distri- 
bution of Vd is Gaussian, with zero mean and variance 



(V d (r,z)V d (r',z')) = A(r-r',z-z') 



(2) 



We next consider the canonical partition function of this 
system (here T is temperature, and we use units such 



that Boltzmann's constant ks — 1), 



N 



(3) 



and average over the disorder by introducing p replicas 
of the above system and making use of the well-known 
replica trick: 



In Z = lim 



Zp- 1 



p~*0 p 



(4) 



upon which we obtain the following, disorder-averaged 
Hamiltonian: 

p N r i a a 2 

s = EE dz l {K (^) + £ nw-wm 

a=l i=l J 

P N 

- ^E E jdzdz'AtfW-r^z-z'),^) 



a,b=l i,j=l ' 



where the superscripts (a, &,...) label replicas. In what 
follows, it will prove useful to write the flux line position 
at height z, r^(z), as the sum 



Ti(z) = r 0l + Uj(z) 



(6) 



where roi = j;Jdz r^(z) is the cm. position of the ith 
flux line, while Ui(z) is the displacement of the ith flux 
line at height z with respect to r$i, and has the following 
decomposition into Rouse modes, [HI 



Ui(z) = Ui(g„)e l9 " z , 



with the Fourier coefficients: 



u;(<7„) = i Jdz u 4 (z)e tq " 



(7) 



(8) 



In keeping with the spirit of the calculation carried 
out in ref. |9|, in this section we shall perform a sim- 
ple perturbative analysis of the physics encoded in the 
Hamiltonian (JSJ and expand H to quadratic order in the 
displacement field. The pure (disorder free) part of H 
yields: 



ff(0) + jj(1) 

■)')) i.r p. 1 - 1 - 1 'rrti.-r-t 



where 



(9) 



(10) 



a=l ijtj 



is the Hamiltonian of a system of perfectly straight flux 
lines interacting through the potential V(r), and where 



p N 



pure 



K(-A +H%u1 a {z)uUz) 



V dz 



a=l i=l J 
P N 

EEE J dz-^utJz)u^(z) (11) 



»=1 1=1 J&i) ' 
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represents the internal modes contribution to the "pure" 
part. In the al 
are given by: 



part. In the above equation, the coefficients fj,^ and n^p 



N 



(i) 
V at 3 = 



u (ij) 



d a d p V(r a 0l ~r a 0j ) 



(12) 



(13) 



In a similar fashion, a Taylor expansion of the disor- 
der part of H to quadratic order in the displacements 
{uf} gives a decomposition similar to the one in Eq. 
namely 



Hdis = H 



(0) 

dis 



H 



(1) 



(14) 



Here, H^ s is given by: 



£r(0) _ 
11 dis ~ 



(15) 



a, b ij 



and is the disorder part of the disorder-averaged Hamil- 
tonian of a system of perfectly straight flux lines in 
a Gaussian random potential with variance A(r) = 

dzA(r, z). On the other hand, S^ g is given by 



H 



(i) 



p N 



~^££ / "dz f dzr\[uUz) [u? (z)-ut (z')]d a dpA(O,z-z') + 

a=l i=l J 

^T,T, [ dz [ dz 'l [«?«(*) - U U Z ')\ i u U z ) - u W)\ d a d A(r a Ql - r a oj ,z z') 

o=l i^j J 

^ £ £ J ^ dZ J ^ dZ '\ - ^aW] [«^W - "S/sOO] 9a^A(r^ - 4, Z - Z') 



a=£b i,j 



(16) 



It is not difficult to see that the first term on the right 
hand side of the above equation represents same replica, 
single-line contributions to the Hamiltonian of the in- 
ternal modes of vortices, while the second and third 
terms represent contributions to HjP coming from same 
replica, different flux lines and from different replicas, 
respectively. 

Collecting all terms, it follows that H can be written 
in the form 



- fH°) 



H = H 



H 



(i) 



where 



u-(0) = &(0) , H v 



eHI) _ Ml) , II 
11 — n pure t n d 



(0) 
lis 

(i) 



(18a) 
(18b) 



In the spirit of ref. Q , we shall derive an effective Hamil- 
tonian for the internal modes of the flux lines by averag- 
ing over the center of mass positions {r^}: 



(H 



(1) 



(19) 



where the average is carried out with statistical weight 
exp(-ff(°)/T)/Z (with Z = Tiie- m ° } / T )). The pure 
part of H u has already been evaluated in ref. Q, with 



the result: 



jjpur 



P N 

£{£/* 

o=l i=l J 
N 

EE/^ 

*=1 3&i) 



1 r /du?\2 
— K 

2 L 

1 /x 



2 TV- 



T <(z). Uj Q (z)}, (20) 



(17) where the "mass" coefficient /i is given by 



(21) 



and where 



,9o(r - O = 4 E E ^ r - r o^ r ' ~ r ^)>° (22) 

9 i=lj#i) 



is the pair distribution function of the two-dimensional 
liquid formed by the centers of mass of flux lines belong- 
ing to the same replica in the vortex liquid. In a similar 
fashion, we show in Appendix lA"lthat H^ ls = (-ff^)o can 
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be written in the form: 



p N 



H u s = ££/ dz[hK(d z u?f+%{u?(z)f 



£ £ / dz J dz - z ') • - 

££ j dz j dz ' 5^ b) (z-z>)u?(z)-u%z>) 



a^b i,j 



(23) 



Here the long- wavelength disorder contribution 5K to the 
tilt modulus of the flux lines is given by 



SK = - 



1 f°° 
— j_ dzz 2 ViA(r,z) 



r=0 



(24) 



and the "mass" coefficients 5/j, are given by 



fdv ^ 50jO6 (r)ViA(r), 

- 1 ^ 6=1 



= (jy _i )dj _ r / dr 5o(r)ViA(r,z) 



(25a) 



Nd i T 



/" dr ff0 , a#b (r)ViA(r,z)^ 



(25b) 



(25c) 



We immediately note that the two-body coefficients 



5(iij(z) and 5fj,\j°'(z) vanish in the thermodynamic 
(N — > oo) limit. In equation l|25al) above, the sum 
J2b=i 9o,ab(r) denotes the quantity: 

p v 
£5o,ah(r) = g {r) + ^ SWb(r) , (26) 
6(/o)=l 



>6) 



6=1 



where 



iV JV 



ff0 ,a^6(r - r') = 1 ££<£(r - r Q ^(r' - <,))o (27) 
P »=i j=i 

is the pair distribution function of the cm. mode of flux 
lines from different replicas. Now, in Eq. I|25a|l . and 
under the assumption of replica symmetry, all pair dis- 
tribution functions go,a^&(r) are equal to the same func- 
tion <?o( r )> and we may replace the sum X)&^a9o.a#6 by 
(p — I) go, with Eq. (|25a|l becoming: 



Sfj, 



P 



dr ffo(r) + (p-l)ffo(r) ViA(r), 



which, in the limit p — > reduces to 
P 



d ± T 



dv 50(r)-fl„(r) VlA(r). (28) 



In the replica-symmetric ground state considered here, 
the diagonal and off-diagonal (in replica space) pair dis- 
tribution functions should be equal, 



.90 (r) =§o(t) , 



(29) 



and we therefore obtain that the correction 5[i identically 
vanishes, which shows that the bare mass p, generated 
by interactions between flux lines is unrenormalized by 
disorder. 

We now can write the following expression for the ef- 
fective Hamiltonian H u — HP ure + H^ ls of the internal 
fluctuations of flux lines in a vortex liquid: 



p N 



a—l i—1 



1 r 



ff «^££ dz o K R (d z ut) 2 + f i(ut(z)Y , (30) 



where we discarded the terms proportional to 1 /N (which 
vanish in the thermodynamic limit, see Eqs. I|25b|) - 
(|25c|0 . and where the renormalized tilt modulus K is 
given by 



K R = K + 6K , 



(31) 



with SK given by Eq. (|24|l . At this stage, an explicit 
expression for the disorder correlator A(r, z) is called for. 
For point disorder such as oxygen vacancies in HTSCs, 
we shall take A(r, z) = A exp(-(r 2 + z 2 )/2£ 2 ), (with 
the understanding that the correlation length £ is much 
smaller that the average intervortex distance a = 1/y/p), 
upon which we obtain the following expression for the 
effective tilt modulus: 



K, 



K 



2ttAo£ 



(32) 



This expression shows that a flux line liquid in presence 
of a weak pinning potential is equivalent to an unpinned 
liquid but with a higher tilt modulus, i.e. that flux lines 
are stiffened by weak point disorder, which is what one 
would expect based on physical intuition. 

For disorder that is correlated along the direction of 
flux lines, e.g. columnar pins, |l9l | the variance A(r, z) 
is z-independent. Upon using for this case a disorder 
correlator of the form A(r, z) = Ao exp(— r 2 /2^ 2 ), we 
obtain 



K 



L 3 A 
3d ± T£ 2 ' 



(33) 



which shows that the tilt modulus is much strongly renor- 
malized by correlated disorder than it is by ordinary 
point disorder, in agreement with the predictions of refs. 
|20l l3l| . The divergence of the right hand side of the 
above equation, however, signals the breakdown of per- 
turbation theory for correlated disorder, which is best 
treated with other, nonperturbative methods [2(| that 
are better suited to strong pinning situations. 
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III. VARIATIONAL APPROACH 



Variation of the trial free energy 



We now generalize the analysis of the previous Section 
to take into account the effect of possible large fluctua- 
tions of flux line trajectories, by using the replica Gaus- 
sian variational approach for elastic manifolds, [l^ For 
compactness, we shall only give the salient features of the 
calculation, and refer the reader interested in more de- 
tails to ref. ^3 j where a similar calculation was done for 
the pure case, the generalization to the disordered case 
being straightforward. 

We start by introducing the following variational 
Hamiltonian: 



H v =H [{T Qi }]+H 1 [{u O i}] 



where H and Hi are trial Hamiltonians for the cm. and 
internal modes, respectively, and are to be determined 
variationally. Although one can, in principle, use a very 
general trial Hamiltonian for the internal modes of flux 
lines of the form 

Hi = EE l dzdz ' P~\<ln)]ta^Mu%{-q n ), 
a,b i,j 

(35) 

the insight we gained from the perturbative solution sug- 
gests the following, simplified form: 



a.i n^O 



(36) 



-T\nZ x + (H-H v ) x , 



(37) 



where Z\ = Tr(exp(— Hi/T)) and (• • *)i denotes averag- 
ing with statistical weight exp(— H\jT")jZ\, with respect 
to the cm. Hamiltonian Hq leads to the result that the 
optimal choice for Hq is given by (we henceforth use the 

shorthand notation J q = J J^ctj; ): 



(34) H = § E { E d ±P T [ G o - G-^qn)] G(q n ) 



i,j n^O 



E / y(q)e iq - (r o-- r o a J ) e -^r^ (0) 

- E / / dz A(q,z)e lq -^~Oe"^^ W 



T 



a,b J 1 



}• 

(38) 



where we defined 4>ij{z) = ([u^(z) — Uj(0)] 2 ). Fur- 
ther variation of the resulting free energy F v = 
—Th\Zx\ H = fj o with respect to G(q n ) leads to the fol- 
lowing result: 



G-\q n ) 



q 2 V(q)g (q)e 



-9 2 G(0) 



d i T 



dz q 2 A(q, z)[l - cos{q n z)]e- q2 ' l > {z V 2dl 



+ '- I dzq z A(q,z) [g (q) + (p- l)g (q))e 



-q 2 ,P(z)/2d 1 



(39) 



r 



where now <fr(z) denotes the relative displacement of in- 
ternal modes within the same flux-line, 4>(z) = ([ui(z) — 
Uj(0)] 2 ), and where the tilde indicates that G~ 1 (q n ) has 
been averaged over all possible configurations of the cm. 
positions jr", | which are compatible with a liquid struc- 
ture. 0, Il0| Under the assumption of replica symmetry 
(which implies that go(q) = go(q)), the last term vanishes 
again in the limit p — ► 0, and we obtain: 



^(qn) 



ql 



?^(q)3o(q)e-' 2( " 2)/dl + 



+ 



d\T 



dzq 2 A(q, z)[l - cos(q n z)}e 



modulus of flux lines in the long wavelength limit: 
1 



K 



R 



K 



d±T 



dz q 2 z 2 A(q,z)e~ q tW/ 2 ^ .(41) 



We thus see that the inverse propagator for the elastic 
distortions of flux lines in the vortex liquid is given by 
the following, generic form: 



G- 1 (q n ) = L(K R q 2 n + m ) 



(42) 



■£■*(>) M 



(40) 



where the effective "mass" coefficient of internal modes 
fluctuations ixr is identical to the quantity derived in ref. 



Mi?. = I q ^(q)g (q)e" 



(43) 



The last term on the right hand side of the above expres- 
sion leads to the following renormalized value of the tilt 



except that it now depends on Kr (through (it 2 )) and 
hence on the disorder strength Aq. 
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In Eq. Q41J1. the relative displacement of the internal 
modes of a given flux line 



(z) = d ± T^2 G(q n )[l - cos(g„z)] 
d i T 



(44a) 



[l - exp ( - y/n/K R \z\)\ (44b) 



depends on Kr, and we therefore see that Eq. i|41[l is in 
fact a self-consistent equation for the effective tilt mod- 
ulus. For the explicit evaluation of hr(T) and Kr, we 
shall make use of the analytical ansatz of ref. [9j for the 
pair correlation function go(j), which is given by: 



g (r) = 1 - T]e^q>(-ar 2 /a 2 ) 



(45) 



where a is a constant of order unity, and < ?; < 1. The 
numerical constant rj quantifies the degree of correlation 
between cm. positions of flux lines. It is close to unity 
when flux lines are strongly anti-correlated due to the 
repulsive interactions between their surrounding super- 
currents, and close to zero in situations where there is 
considerable cutting and crossing of flux lines. Using the 
above ansatz for ga(r), we obtain the following expression 
for the "mass" /j,r of the internal modes as a function of 

t, E3 



Mo 




uT 



Aa 2 ^K R ^ ) \4a 2 y/K R fi J 



aT 



(46) 

where /io = Hr{T = 0) = 2riirpEg/d \ . Since fiR is a 
monotonically decreasing function |10| of the parameter 
v = (aT/4a 2 Kr/iq ) , we arrive at the important con- 
clusion that pinning disorder (which increases the value 
of the tilt modulus from its bare value K to the effective 
value Kr > K) increases the value of /i, thereby reduc- 
ing even further (than the sole increase in K) the thermal 
wandering of flux lines. This effect, which did not appear 
in the elementary treatment of Sec. |nl is expected to be 
rather weak for the weak disorder considered in this work, 
but may nevertheless reveal itself in actual experiments. 

Going back to the effective tilt modulus of the pinned 
flux line liquid, we see that due to the highly nonlin- 
ear character of the self-consistency equation H41() . it is 
not possible to solve this equation for Kr and obtain a 
general expression for the effective tilt modulus in closed 
analytic form. For weak disorder, however, such that 
Ao <C KT/£, one can evaluate the second term in Eq. 
(|41|l perturbatively in Ao, by using for <fi(z) its expression 
in the absence of disorder, Eq. (|44b[l . with K instead of 
Kr. In the limit of small correlation lenght £, we obtain 



Kr = K 1 



v^A xe 3 

2V2d ± T 2 



(47) 



Comparison with Eq. (|32fl shows that the correction to 
the bare tilt modulus K goes to zero at high temperature 
more rapidly than in the simple Taylor result of Eq. Il-Sl'l) 



(1/T 2 as opposed to 1/T). This discrepancy is due to 
the fact that the Taylor expansion of Sec. [n] does not 
take proper account of thermal fluctuations of the inter- 
nal modes of flux lines, as opposed to the self-consistent 
approach of this Section which in fact can be shown [l5j 
to become exact in the limit d j_ — > oo. 



IV. ACTION FORMULATION OF THE 
DYNAMICS OF FLUX-LINE LIQUIDS 

We now turn our attention to the dynamics of vor- 
tex liquids. The general dynamical behavior of flux line 
assemblies (solids and liquids) in a random pinning envi- 
ronment has attracted a lot of attention in recent years 
due on the one hand to the considerable technological im- 
plications of understanding the flow of vortices in HTSCs, 
and on the other hand because of the fundamental theo- 
retical questions and variety of physical regimes displayed 
by these systems. Most of the above mentioned attention 
has focused on the dynamics of flux solids, with ques- 
tions about the degree of crystalline and temporal order 
in the driven regime, and glassiness in the absence of 
external drive, at the forefront of theoretical issues that 
have been addressed. Although there have been a num- 
ber of studies of the dynamics of disordered flux liquids, 
these studies were either done within hydrodynamic ap- 
proaches which, by definition (since they involve a coarse 
graining procedure over many vortices) are unable to give 
information about actual flux-line trajectories inside the 
superconducting sample, or used qualitative arguments 
to separate length and time scales in the plastic regime 
near the melting point. 14] While the latter approach is 
very valuable in that it helps draw a good qualitative pic- 
ture of the physics of driven flux liquids, it assumes that 
the vortex liquid is very viscous and hence only applies 
very close to the melting point. As we mentioned in the 
Introduction, here our goal is to go beyond these previ- 
ous treatments, and establish a general framework for a 
systematic perturbative study of driven flux-line liquids 
in presence of weak disorder. 

We shall assume that the motion of flux lines in the 
liquid state in presence of an external driving force F is 
governed by the following, overdamped Langevin equa- 
tion: 



5H 



8ri(z,t) 



F + Q{z,t). 



(48) 



For simplicity, the distribution of the thermal noise 
Ci{z, t) will be taken to be Gaussian, with zero mean and 
correlations (we use units such that fc^ = 1): 

(C ia {z, t)C jP (z', 0) = 2 7 T%<5 Q/3 S(z - z')5(t - t'). (49) 

In the above equations, the parameter 7 is the micro- 
scopic friction coefficient characteristic of the interaction 
of the system with the degrees of freedom of the sur- 
rounding heat bath. In our case of flux lines in a type II 
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superconductor, 7 describes the drag on a flux line due 
to the interactions of the electrons in the normal vor- 
tex cores with the underlying solid, and is given by the 
Bardeen-Stephen expression: [lT], |2l| 

7=8^-^, (50) 

where p is the average density of flux lines, £ a b is the 
superconducting coherence length in the (ab) planes, and 
(T,i is the normal-state conductivity. 

The dynamics represented by Eq. (|48|) is best stud- 
ied u sing the action formulation of Martin, Siggia and 
Rose |H |H (MSR), whereby disorder-averaged 

observables are obtained from the following, disorder- 
averaged generating functional: [8j 

N 

Z = J\{[dv l {z,t)][dv l {z 1 t)\e~ A . (51) 

i— 1 

Here, the MSR dynamical "action" can be written in the 
form: 

A -^/ree ~\~ Aint ~t~ -Adis j (52) 

where the "free" part 

Afree = E / dz dt [- (2^)*? (z , t) 

+ iii(z,t)- [jd t ri(z,t) - Kd 2 z Ti{z,t)]} (53) 

corresponds to an "ideal gas" of non-interacting flux 
lines, and 

■Aint=^2[ dzdtir l (z,t)-VV(r l (z,t)-r J (z,t)) (54) 

is the part of the action describing the interactions be- 
tween vortices. The last term in Eq. i|52|) decribes the 
pinning of flux lines by the underlying disorder potential, 
and is given by: 

Adis = — — 5^ fdzdt I dz 'dt' ' r na (z,t)r m/3 (z' ',t') 
n,m 

x d a d p A{r n {z,t)-v m {z',t')). (55) 

(In the above equations, and in what follows, Y] i < ■ 

stands for the double summation J^Li Sj(^j))- Our 
main goal in the next two Sections will be to find a way 
to calculate expectation values of dynamical observables 
in the liquid phase (like for example the dynamic struc- 
ture factor of the flux line liquid) by integrating directly 
over the conformation variables {?i(z, t), Vi(z, t)} using 
the MSR generating functional of Eq. i|51|) . and not by 
integrating over the density operators, as is done in the 
(static) boson analogy and other hydrodynamic ap- 
proaches. 8] This means that we have to avoid writing 



the dynamical action A in terms of the density opera- 
tor p(r, z; t) — J2iLi <H r — v i{ z i *))) an d instead keep the 
conformation variables {f i(z, t), rj(z, t)} as the true and 
only dynamical variables in the problem. Our strategy 
will be very similar to the strategy adopted in the pre- 
vious two Sections, which consists in separating out the 
center of mass from the internal modes of flux lines, and 
trying to find a decoupled approximation to the dynamic 
action A 

A = A {0) + A {u) , (56) 

such that A^ and A^ depend only on the cm. and on 
the internal modes, respectively. Since our ability to per- 
form functional integrations is limited to Gaussians, for a 
decomposition of the form (|56[) to be useful at all we will 
need to write A^ as a bilinear form in the displacement 
fields {u.i(z, t), Ui(z, t)} fields. Such a decomposition will 
allow us to evaluate averages of the form: 

(p{v,z;t)p(r',z';t')) = ^ (S(r - r 0n {t) - u n (z, t)) 

x S(r' -r 0m (t')-u m (z',t'))) 

(57) 

by integrating over the cm. modes {roi(i), foi(i)} and 
internal conformation variables {ui(z, t), u^(z, t)} of flux 
lines, which are the true dynamical variables of the 
system, rather than the averaged density p(r, z;<) = 
(p(r,z;t)} which is a coarse-grained density with no de- 
tailed information on these conformation variables. We 
shall first start by considering the case of a pure (disor- 
der free) flux liquid at equilibrium, which will serve as 
a starting point to our treatment of the disordered case 
to be considered in Sec. IVII and needs therefore to be as 
accurately specified and perfectly understood as possible. 

V. EQUILIBRIUM DYNAMICS OF PURE FLUX 
LIQUIDS 

In the free part Af re e of the action, Eq. I|53[l . we 
rewrite the flux line trajectories {ri(z,i)} and conjugate 
fileds {fj(z,£)} in the form: 

Vi{z,t) = r 0i (t) + u;(z,t) , (58a) 
ii(z,t) = r i{t)+Ui(z,t). (58b) 

Using the fact that J dz Ui(z,t) — J dz Ui(z,i) = 0, we 
easily obtain that the free part of the action Af ree can 
be written as the sum: 

Afree = ^/ree + ^free ' (59) 

where 

A free = J2 / dt {-{2^)^) + i^t) L^V^t)) 
i=l •* 

(60) 
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depends only on cm. variables, while 



with 



J dzdt{\{2 1 T)u 2 i {z,t) + 

i—i 

iiii(z,t) [yd t Ui(z,t) - Kd 2 z Ui(z, t)] } 



is the free action for the internal modes of the flux lines. 
From Eq. (|6U[) . it is not difficult to see that the cm. 
mode of flux lines is characterized by a friction coefficient 
70 = Lj, and hence that the diffusion constant Do of a 
free flux line is inversely proportional to the thickness L 
of the sample, as already pointed out in refs. |sl l26|: 



Do 



1 

17 



In the interaction part of the action Aint , we expand the 
interaction potential V(ri(z,t) — Tj(z,t)j to linear order 
in the displacement field: 

V(r t (z,t) - rj{z,t)) = V(r 0l (t) - r 0j (t)) + 

+ [ui(z,t) - Uj(z,t)} ■ VV(r 0i {t) ~ r 0j {t)). 

The interaction part Ai n t can then be written in the form: 
Aint — A int + A\ n \ , (62) 



-4n* = E/ dtiLr i(t)-W(r 0i (t)-r 0j (t)), (63a) 

•Ant = E J dzdtiui a (z,t)[uip(z,t) — Ujp(z,t)]x 

x d a dpV(r 0i (t)-r 0j (t)). (63b) 

Combining Eqs. 1(59(1 and 1(62(1 . we see that we can already 
write the total action A in the form A = A^ + A^\ 
where 



(61) A® 



A 



(0) 



free ' int 



A 



(0) 



N r 1 

E / rft {-(2£ 7 T)f^) + »«(*) • [i7^r i(t) 



+ LWV(r 0i (t)-r 0j {t)) 



(64) 



depends exclusively on cm. variables and can be thought 
of as the dynamical MSR action of a liquid of hard rods 
of length L interacting through the potential Vo(r) = 
LV{v). The effective action A^ is on the other hand 
given by: 



N 

A [1) = J2j dzdt{~(2'yT)u i (z,t)+iUi(z,t)- [jd t vu(z,t) - K d*va(z, 



t) + 



i=i 
iu 



a (z,t)\( ^2 d a dpV (r 0l (t) -r ok (t))^j 5 tj - ^ d a dpV (r 0i (t) -r 0j (t)) Ujp{z,t)\ 



and describes the internal fluctuations of the flux lines. 
As it stands, however, A^ still contains cm. dynamical 
variables. In order to obtain an effective action which 
depends only on the internal modes, we need to take the 
average of A^ over all configurations of the cm. coordi- 
nates {r oi (t)} which are compatible with a liquid struc- 
ture. In the spirit of a cumulant expansion, |9j we shall 
write A^ = (A^)o, where the average here is taken 
with statistical weig ht exp(-^°)), A w being the dy- 
namical action for the cm. mode, Eq. 164(1 '). Performing 
the above Gaussian average, we obtain: 



E / &d^(2 7 T) Ul 2 (M) 
i=i J 

+ iu ia (z,t) U^dt- Kdl)5 a p + ^Auip{z,t) 

+ E ^S^/^)]}' ( 66 ) 



(^i(r-r oi (f(r'-r oj (t))) o; 



(65) 



where we defined: 9] 

= (E 9 « 5 ^( r «(*)- r »i(*))) > ( 67a ) 

^? = -(d a dpV(r 0i (t)-r 0j (t))) Q . (67b) 
It is easy to see that 

fj.% = J dvdv' d a d{iV(v-v') x 



p 2 J drdr' d a d p V(r - r') g (r - r') , (68) 



where 



,9o 



( r - r ') = 3 < E <*( r - *<*®W ~ *0j®))a (69) 
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is the (equal time) equilibrium pair distribution function 
of the cm. mode of the flux line liquid. Using the rota- 
tional symmetry of both go(r) and V(r), we obtain that 

/i = j- J dr g (r)V 2 V(r) , (70) 
which is the result l|21|) of Sec. [HI and similarly that |j| 

llSa/3 



N-l 



(71) 



The above results lead to the following expression for the 
effective action of the internal modes of flux lines: 



i—l 
N 

- ^2 iui a {z, t) (^yd t - Kd 2 z + N _ 1 A* )<■> 



N- 



-] Uj{z,t) | 



(72) 



The above effective action can be written in Fourier space 
in the form: 

A (u) = Y^Yl j [^ fl i(9n> W )fy(9n,w) • Uj(-q n ,-Uj) 



+ iUi(q n ,w)Tij(q n ,uj) ■ Uj{-q n ,-u) 



(73) 



where the shorthand f stands for f°° and where 



the dynamical kernels r^cfo, w) and Tij(q n ,u>) are given 
by: 



Tij(q n ,Lu) = 2^/TLSij, 
Tij(q n ,uj) = L ( - £70; + ATg£ + 



(74a) 



iV 



TV- l 



iV- 1 

In the thermodynamic limit N — > 
Tij(q n ,ui) reduces to the diagonal form: 



(74b) 

00, the kernel 



(75) 



We hence obtain in our perturbative approach that the 
internal modes of different flux lines are effectively de- 
coupled: in the thermodynamic N — > 00 limit, the effect 
of the interactions between vortices on their internal fluc- 
tuations is entirely encoded in the /1 term, which acts as 
a quadratic confining potential (~ \l^u 2 ) for the internal 
modes of individual flux lines (in total agreement with 
the findings of the static approach of ref. 0|). 

Knowledge of the propagator Tij(q ni Lu) allows us to 
find the reponse and correlation functions, Rij (q n , ui) and 



Cij{q n ,uj) respectively, which we define as follows: 



Rij{z- z',t-t') 



C ij (z-z , ,t-t') 



6(u ia (z,t)) 
5Q a (z>,t>) ' 

— (m (z,t) -iuj (z',t')), (76a) 
dj_ 

-^{u l {z 1 t)-u. J (z'X)). (76b) 

dj_ 



Whithin the mean-field approach of this Section, and in 
the thermodynamic limit N — * 00, the above functions 
are both diagonal in the vortex labels i, j. If we denote 
by T(z,t), R(z,t) and C(z,t) the diagonal parts of the 
vertex, response and correlation functions respectively, in 
such a way that [z, t) = T(z, t)Sij , Rij [z — z' , t — t') — 
R(z,t)Sij and G'g(z — z\t — t') = C(z,t)Sij, then one can 
easily verify that [H El IH 



R{q n ,u) 



(77a) 
(77b) 



In the thermodynamic (N — > 00) limit, we obtain from 
Eqs. lj?5)> and Ij77a|> that the response function R(q n ,u) 
is given by: 



R(q n ,u) 



L(-ijuj + Kql + n) 



q n ^0- (78) 



Performing a partial Fourier transform back to the vari- 
able t, we obtain: 



R(q n ,t) 



(79) 



where 9 is Heaviside's unit step function. On the other 
hand, from Eq. I|77b|) . we readily obtain for the correla- 
tion function C{q n ,t) the following expression: 



C(?„,t) 



T 



L(Kql + f i) 



-(p+Kql)\t\h 



q n ^ 0. 



( 8 °) 

It is easy to verify that the fluctuation-dissipation rela- 
tion 



6(t)d t C(q n ,t) = -TR(q n , t) , q n ^ 



(81) 



holds for the internal modes of flux lines, which indicates 
that these modes will eventually reach thermal equilib- 
rium at long enough times. We however should empha- 
size that the above expressions of the response and cor- 
relation functions are only valid for the internal modes of 
the flux lines. The cm. mode of vortices, as described by 
the effective action (|64f) . is still diffusive, although with 
a diffusion constant D which we expect to be reduced by 
the interactions to a value that is smaller [2]} than the 
bare diffusion constant of free, noninteracting vortices 
Da = 1/(7 L) of Eq. (|61|l . Indeed, from ref. [27|, one can 
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estimate the value of the interacting diffusion constant 
D in the absence of disorder to be of order: 



1 + 2^ 2 <?o(2a) ' v ' 

which is smaller than the bare diffusion constant Dq for 
all values of applied magnetic fields smaller than the up- 
per critical field H C 2- 

Having derived the equilbrium dynamics of pure flux 
liquids, we now turn our attention to the more general 
case of driven vortex liquids in the presence of a ran- 
dom pinning potential. It should be pointed out at this 
stage that the "mass" coefficient /i given in Eq. H70|) , can 
be generalized to take into account large vortex distor- 
tions, as was done for the static case in ref. [lOj. This 
generalization is performed within a dynamic Hartree ap- 
proximation in Appendix B. 

VI. DYNAMICS OF DRIVEN FLUX-LINE 
LIQUIDS: PERTURBATION THEORY 

A. Coarse-grained effective action for flux line 
dynamics 

Following refs. 0, H^ . Il3| , we decompose the internal 
modes of flux lines into short- and long-wavelength parts, 

Ui{z,t) = u<(M) + u>(z,i), (83) 



where (we here for convenience adopt a continuous nota- 
tion for the q n summations): 



u<(z,t) = f u t (q z ,t)e^ z , (84a) 
u>(M) - f Ui(q t ,t)e^". (84b) 

JA < <q<A 



In the above equations, A and A< are high and low mo- 
mentum cut-offs, respectively. The ultraviolet cutoff A is 
given in terms of the superconducting coherence length 
£ c along the direction of the flux lines by A = 7r/£ c . In- 
serting the above decomposition, Eq. (|83[) . and a simi- 
lar decomposition for the response field Uj(z, t), into the 
dynamical action, we find after integrating out the short 
wavelength modes {uf , uf } that the long- wavelength ef- 
fective action, to first order perturbation theory, is given 
by: 



Aeff = Apure + {A d is [H* + U> , U< + U>] ) > , (85) 



where 



(A dis ) > = ^^[dzdtfdz , de?*{z,t)^(z,t;z',1?)T%(z , J^ (86) 

i,j i 

where we defined 

STf(z 7 t;z',t') = /i ?a j^A(q 1 z-z')e^> T « < ^- , J < t , '^e-i to w*Jf(^' t '), (87a) 
Sr ia (z,f,z',t') = I dz'dt'Rijiz^iz'^') f g a <z 2 A(q,z-z')e iq ' [rf(z ' tKrj<(z '^' )l e"^ Q9 ' 5 ^' 3(z '' ;z '''' ) (87b) 

with (j)"^ (z,t; z',t') — ([uia(z,t) — Ui a (z', t')][v,jp(z, t) — Ujp(z' ,t')]) . In perturbation theory, valid at large driving 
forces, it is convenient to take the limit A< — > 0, and use for the response and correlation functions the expressions 
(|79fl - l|80ll . which are spatially invariant and diagonal in vortex indices, in which case expressions (|87all - (|87b[) are 
greatly simplified, and become: 

5T?Hz,t;z>J) = / l(Z ^A(q,z~z')e iq -[ r ^ 2 ^- r , < ( z ''*')]e-^T* (z - z ''*" t ' ) , (88a) 



ST?(z,t;z',t') = ^ j dz'dt' R(z- z';t-t') j iq a q 2 A(q, * ~ z 1 )^^^-^'^-^^-*'*-^ 



I 

In the above expression, $(z, t) denotes the quantity where (f>o(t) = ([r (i) — r (0)] 2 )o is the relative displace- 

ment of the cm. mode, while 

*(>,*) =0o(i) + 0OM) (89) 

<j>{z,t) = ([ Ul (z,i)- Ul (0,0)] 2 ), 

2Td ± ^ 1 - cosfe^e-^+M) 1*1/7 
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is the relative displacement of the internal mode u of a 
given flux line in the vortex liquid. 

Using Eqs. (|S51) . (JHUJl and (|88b|> . effective equations of 
motion for the cm. and internal modes of flux lines can 
be derived in a standard way. [Til IT2I Il3| For the internal 
modes, we obtain 

l a pd t u ij3 {z, t) = (-(J,8 a p + K af3 d 2 z )u ifj + 

- ^aPfdziiipdzUij + Q{z,t) , (91) 

where now, in addition to the usual (but renormalized) 
elastic tension term, new non-linear Kardar-Parisi-Zhang 
(KPZ) terms have appeared. The cm. mode on the other 
hand obeys the following equation of motion: 



lapdtroipit) = F a - Ffr,a 



(92) 



where the friction force Fyy arises from the additional 
drag experienced by the vortex liquid as a result of the 
presence of the random pinning potential. In the follow- 
ing subsection, we outline the main steps of the coarse 
graining procedure which leads to the above equation of 
motion, before moving on in Sec. IVlll to calculating the 
dynamic structure factor of the flux line liquid in pres- 
ence of disorder. 



2. Friction coefficient 

The disorder-correction to the friction coefficient is ex- 
tracted from the effective action as follows. In the ex- 
pression (|88b|) of the effective kernel ST ia , we make use 
of the expansion 

e iq-[u 1 («,t)-u,(*',t')] ~ l + i qa [u ia (z,t)-U ja (z',t')} 

(96) 

We further shall assume that the disorder is weak, so that 
the internal displacements of flux lines vary slowly on the 
scale of the equilibrium kernel T(z, i), i.e. on the scale of 
y^n/K. In this case, one can make use of the following 
gradient expansion: 

Ui a (z,t) - Uj a (z',t') ~ {t - t')dtUi a (z,t) + 

+ (z- z')d z u la - l -(z- z'fd 2 z u la . (97) 



Using both expansions and lj9"7|l into Eq. (|88 b|) . we 

obtain that the effective action A e ff contains a term of 
the form: 

Aeff = J dzdt iff a (z,t)[S-f al3 d t u if3 (z,t)} , (98) 



B. Derivation of renormalized quantities 

1. Friction force 

The first order correction to the friction force is ex- 
tracted from the ~ 0(?oi) terms of the effective action 
ofEq. (JHH): 

Aeff [roi] = -J^Yl j dtLi 0ia{t) J dz J fc? 2 A(q, z) 

x i?(z,i)e 4qvt ~^* (M) . (93) 
This term is of the form: 

A eff [r 0i ] = Y:JdtiLr 0ia m r , a , (94) 
with the effective friction force: 

F fr , a = Jdzdt I ^-A(q, z)R{z, t) sin(q • vt)e~^I Hz ' t} . 

This expression is identical to previously derived expres- 
sions |l4| for the friction force on individual flux lines in 
presence of disorder, except that here the response and 
correlation functions to be used for an explicit evaluation 
of -F/r,a should be ones that are relevant to a flux liquid, 
e.g. Eqs. (O and JSSJl respectively. 



and hence that the effective friction coefficient jap in Eq. 
lHIH is given by: 



7q/3 = 7<W + , 



(99) 



with 



S laf3 = Jdzdt J ^-A(q,z)^(z,i) e iqvt -^r* (2 '* ) . 

(100) 

Note that in the absence of drive (v = 0), 7 Q/ 3 is isotropic, 
7a/3 = (7 + 6j)6 a p, where now: 



6-y = -p- I dzdt I q li A(q 1 z)tR{z,t)e" 



3. Elastic dispersion 



(101) 



It also follows from Eqs. (|96[> - (|97fl that the elastic co- 
efficients appearing in Eq. (|91|l can be written in the 
form 



K a p = K5 a f3 + SK a /3 , 



(102) 



where the disorder-dependent correction 8K a p is given 
by: 



5K af} = J dzdt J q a q l3 q 2 z 2 A((i,z)R(z,t) x 



(103) 
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This equation can be rewritten, using the fluctuation- 
dissipation theorem, Eq. (|81|) . in the form 

SK af:} = - J dz J q a q/3Z 2 A(q,z) x 

x / dt e^-^Z^^dte-™!^ (104) 
Jo 

In the static (v = 0) limit, the above correction to the 
tilt modulus becomes isotropic, bK a p — SKSap, with: 

5K = ^- f JdzJ g 2 z 2 A(q,z) e ^ $ ^+0(l/L), 

(105) 



which is the result (|41|l that we obtained in Sec. IIIII 
within a static replica approach. 

The Taylor expansion l|97|l also yields convective terms 
of the form Xapcl z Uip{z, t) on the right hand side of Eq. 
(|9ip. It is easy to see, however, that the coefficients of 
these convective terms, 



Xa/3 



— J dzdt J q a qpq 2 zA(q, z)R(z,t) 



2d 



x e *»- vt e-*r* ( *' t) , 



(106) 



vanish identically by virtue of the fact that the integrand 
is odd in the integration variable z (provided that the 
disorder correlator is even in z, i.e. A(r, —z) = A(r, z)). 



4- The KPZ nonlinearity 

The effective action (|85|l contains an additional contri- 
bution of the form 



A e ff[uu 2 } = Ylj dzdtirf a (z,t) 



X { - -Xayp[d z Us(z,t)][d z Up(z,t)]}, 



(107) 



with 



d± 



J dzdt J q a qpq^q 2 z 2 A(q, z)R(z,t) 



(108) 



This means that disorder induces a KPZ nonlinearity in 
the driven state, much as it does for driven vortex solids. 
In fact, the emergence of non-linear KPZ terms in driven 
flux-line liquids has been predicted a long time ago within 
a macroscopic approach by Hwa |32^ . who found that 
these terms affect the dynamics of the flux liquid on long 
length scales, with the vortices forming a smooth, lami- 
nar phase at small drives, and a rough, turbulent phase 
a large values of the applied force. We shall not study 
the effect of the KPZ terms in any more detail here, and 
instead refer the reader to the above mentioned paper for 
more details on this particular question. 



VII. DYNAMIC STRUCTURE FACTOR OF 
DISORDERED FLUX LINE LIQUIDS 

We are now in a position to calculate the dynamic 
structure factor S(r, z; t) of our flux line liquid. By defi- 
nition: 

S(r,*;i) = (#r,z;t)p(0,0;0)>, (109) 

where (• • • ) now stands for averaging over both cm. and 
internal conformation variables of vortices, and where 
space translational invariance of the flux-liquid has been 
assumed. In what follows, it will be convenient to con- 
sider the partial Fourier transform 

S(q, z- 1) = -4l(p(q, z; t)p(-q, 0; 0)) , (110) 
L ± 

where L± is the size of the system in the plane perpen- 
dicular to flux lines. Using the definition of the density 
operator at time t, 

N 

p(r,z;i) = ^<5(r-r,(z,<)), (111) 

i=l 

we readily obtain that the Fourier transform S^q, z; t) is 
given by 

N N 

S(d, Z ;t)4EE \ . (112) 

L i i=l 3=1 

We now separate the cm. mode from the internal modes 
of the flux lines, and assume that the dynamical action 
A has been written in the decoupled form: 

A = A {0) +A {u) (113) 

where ^4^°^ depends only on the cm. variables {i*oj(t)}, 
while A^ only depends on the internal modes {u^(z, t)}. 
It then follows that the average on the right hand side of 
Eq. (|112f> can be written in the form: 



e -iq-[r i (z,t)-r 3 (z',i')]\ _ /g-iq- [r i(t)-roj (*')] 





5i-, 2 ([u,(z,t)-u,( 2 ',t')] 2 ) 



x e 



(114) 

In the above expressions, (• • • )o and (•••)« denote aver- 
ages with statistical weig hts exp(-^(°)) and exp(-^( u )), 
respectively. Now, in the approaches of Sees. HVI and lVll 
the internal modes of different flux lines are decoupled, 
which implies that, for i =/= j, 

([u i (z,t)-u j (z',t')] 2 )=2(u 2 (z,t)), (115) 

and so the expression (|112|l of 5(q, q z ; t) becomes 

J- i— 1 



^^^e-iq-IroiM-ro^O)]^ g - ^ <« 2 (*,«)> j 

(116) 
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Given that all vortices in the flux liquid are equivalent to 
each other in our mean field approach, we see that the 
first term on the right hand side of Eq. (|116|l involves 
the sum of N identical terms, which we can simply write 
as NF(q,t), with: 29] 

i N 



i=l 



e * a x (117) 

On the other hand, it is easy to verify that 

o 



^^e-^ [ro<(i) - r « (t ' )] } n = L^p 2 5 o(q,t-t'), (118) 



where go(q; t — t') is the Fourier transform of the time de- 
pendent version of the pair distribution function of stan- 
dard liquid state theory, which, in real space, is given by: 



5o (r- 



v 0n (t))S(r'-v Qm (t'))) Q . 

(119) 

Using the fact that N = L d ^ p, we finally obtain 
S(q, z; t) = pe-^I q2 * {t) + p 2 g a (q, t - t') e~t 

(120) 

A principal difficulty with the expression above is to find 
a good approximation for the time-dependent pair corre- 
lation function <?(q, t — t'). The simplest such approxima- 
tion is the so-called "convolution approximation" , which 
was introduced many years ago by Vineyard, |28( and 
which consists in writing for g (q, t) the following ex- 
pression: 



5o(q,0 = go{q)F(q,t) . 



(121) 



Here <7o(q) = go(q, t = 0) is the ordinary (equal time) 
pair distribution function of static liquid state theory, 
and F(q,t) is the quantity defined in Eq. (|117fl . Using 
this approximation into Eq. I|120|) above leads to the 
following result: 



S(q, z; t)=pe-^I Mt) {pgo(q)e-^ {u2 



+ e 



0(2,*) I 
(122) 

where we used Eq. (|117fl to express F(q, t) in terms of 
4>o(t). In the case where the cm. mode of the flux lines 
obeys a simple diffusion law of the form: 



bo(t) = 2d±TD\t\, 



(123) 



we obtain the following expression for the structure factor 
S(q,z;t) 



S(q,z;t) = pe- DT ^ 2 L 50 (q)e 



(124) 



In the following, last Section, we shall discuss some limit- 
ing cases, and compare our theoretical prediction for the 
static structure factor to experimental results. 



VIII. DISCUSSION AND CONCLUSIONS 

We now discuss the meaning and phsical implications 
of our results. We shall start by addressing the nature of 
the unusual massive phonon mode of the internal fluctu- 
ations of flux-lines, which seems to violate translational 
invariance. The first occurrence of this massive mode is in 
Eq. (|llf) . and it is very easy to verify that this equation is 
translationally invariant (even though this is not obvious 
from the way it is written). In Eq. (|30|l . the off-diagonal 
p,ij terms, which are necessary to maintain translational 
invariance, were discarded for simplicity. Keeping these 
off-diagonal terms only leads to corrections of order 1 /N 
to the elastic propagator (as was shown in detail in ref. 
9]) which vanish in the thermodynamic limit. In fact, 
even though a massive term may seem unusual, it is very 
well known that massive phonon modes do appear in or- 
dinary crystal lattices if the lattice has a basis (these are 
the so-called [33| "optical phonons"). In the case of flux 
liquids, the internal modes of the flux lines are the ana- 
logues of atoms belonging to the basis, and the p, terms 
may be thought of as the analogues of "optical phonons" 
in crystals. 

We now turn our attention to the static structure fac- 
tor S(r, z) of the flux line liquid, which is defined as: 

S(r,z) = (p(r,z)p(0,0)) ■ (125) 

In ref. we have shown that for an effective Hamil- 

tonian H u of the general form: 

1 N 

i=l n^O 

the structure factor is given by: 

S(q, z) = p 2 g (q) e^ 2 <« 2 )/^ + ^-dr^W ) (i 27 ) 

where the mean square relative displacement <fr(z) has 
already been evaluated in Eq. i|44b|) , and where the mean 
projected area (u 2 ) of a given flux line is given in terms 
of the elastic propagator G by: 



(u 2 ) =d 1 _TY,G(q n ) 

n=£0 



(128) 



In the present case of a disordered vortex liquid with 
the inverse elastic propagator of Eq. (|42|l . the above 
expression for (u 2 ) yields: 



(129) 



2y/p; R K R 

We thus obtain for the structure factor the following ex- 
pression: 

Tq 2 



S(q,z) 



P <Mqj exp 



) 



pexp • 



Tq 2 



2^/p R K R 



(130) 
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The above result for the structure factor is very different 
from the result obtained by previous authors in the hy- 
drodynamic formulations of refs. IE El IE 01 which is 
given by: 



S(q, q z 



pTg 2 /K 
gf + e 2( q )/ T 2 



A(q,g z ) 



pq 2 /K 



e(q) _ 


\(—Y- 


pV{q)q 2 - 


1/2 


T 


A 2KJ 


K - 





-<?2+e 2 (q)/T2 

where the excitation spectrum has the usual bosonic |3 
form: 



(132) 



In ref. [ToT ] we made a number of comments on the struc- 
ture factor of interacting, but otherwise disorder-free flux 
line liquids derived in the boson mapping approach, Eq. 
(| 1 3 1 11 , and revealed that it had a number of quite disturb- 
ing inconsistencies. We here make the similar observation 
that while the form (|127fl follows from the very general 
assumption that a decomposition of the form l|34|l (with 
an arbitrary elastic propagator G(q n ) in Eq. (|36|l ). can 
be written for the Hamiltonian of the vortex liquid, the 
structure factor of Eq. (|131J) does not correspond to any 
choice of elastic propagator G(q„), and cannot be derived 
from a microscopic approach like ours. In Appendix IUlwe 
show that an expression for the structure factor that is 
similar to Eq. (|131|l (with an additive contribution pro- 
portional to disorder) can be derived for the correlations 
of the density of the center of mass coordinates of flux 
lines. One may therefore speculate that the discrepancy 
between our results and those of refs. 0, 0, 0] are due 
to the fact that in these previous studies the non-diffusive 
internal modes of vortices were treated on equal footing 
with the cm. mode (a diffusive mode for which hydro- 
dynamics is naturally expected to be valid). A more de- 
tailed discussion of the formulation of Gaussian hydro- 
dynamics of continuous media and of the shortcomings 
of previous attempts to formulate the hydrodynamics of 
flux-line liquids can be found in Appendix [Cl 

We now consider some limiting cases. For a "perfect 
gas" of noninteracting flux lines, go(r) — 1, and Eq. 1122|) 
reduces to: 



S(q,z;t) = (27T) 2 p 2 S(q)+pe 



-D T\t\q 2 



(133) 

with the diffusion constant of the noninteracting cm. 
mode Dq = 1/ ("/L) . For noninteracting flux lines, the 
correlation function 4>(z,t) is given by: 



CM) 



2Td 



U n=l y " 



At times shorter than the characteristic Rouse time 
tRouse = jL 2 1 ' K, the sum can be transformed into an 
integral, with the result @: 



2Td 1 
K 



(135) 



where / is the function given by (here T(a, x) is the in- 
complete gamma function [l6l|): 



n JO 

1 

2 " 



dx 



x 
1 



4^ 



r 



1 1 

2' 4u 



(136) 



with the limiting behavior f(u) ~ 1/2 for u — > and 
f(u) ~ \/u/tt for u ^> 1. On the other hand, at long 
enough times, t > tR OUSe , the sum in Eq. l|134|) is rapidly 
converging to the limiting value: 



2Td± 1 
9{t > t Rouse ) = — — > — , 
LK ^ qZ 

LTd± 
12K ' 



(137) 



It thus follows that, at long times (t > t Rouse ), the struc- 
ture factor of noninteracting flux lines can be written in 
the form: 

S(q,z;t) = pe~ D ° T ^i 2 [l + (2^)V(q)]e-^ 

w {2n) 2 p 2 5{q). (138) 

Let us now consider the other limiting case of a liquid of 
infinitely rigid flux lines. If we formally let K —> oo in 
our equations, it is not difficult to verify that (u 2 ) = 
and 4>(z,t) — ([ui(z,t) — Ui(0,0)] 2 ) = 0, and thus Eq. 
(|122fl reduces to the appropriate expression for a liquid 
of point particles, in this case the liquid formed by the 
centers of mass of the different vortices: 



S(q,z;t)=pe-^I q2Mt} [l + pg (q)] 



(139) 



Note that the z dependence has dropped from this last 
equation. Now, since , by definition, 0o(O) = 0, we see 
that the equal time structure factors S(q,z;t — 0) of 
Eqs. I|122|) . I|124|) and l|139|) reduce to the corresponding 
quantities obtained in refs. 0,E1- I n particular, for a 
liquid of rigid flux lines, Eq. I|139|) reduces to the correct 
expression of the static structure factor of standard liquid 
state theory, 



S(q,z) = p[l + pg (q)] 



(140) 



an expression which cannot be reproduced using boson 
mapping and other hydrodynamic methods. 

We now turn our attention to the interacting structure 
factor of Eq. (|122fl . By contrast to the case of non- 
interacting flux-lines treated above, we see here that due 
to the presence of the confining term p,, the relaxation of 
the internal modes is extremely fast: on time scales larger 
than the characteristic time = 7///, the correlation 
function <j>(z,i) of Eq. H9l)fl reaches its limiting value: 



(z,t> tp) 



2Td 



1 



n 

Td ± 



dq 



Kq 2 + p 



(141) 
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which implies that the low frequency behavior of the cor- 
relation function S(q, z; u>) is given by: 



S{q,z;uj) 



2pDTq 2 
• 2 + DT 2 q 4 



1 



•pso(q)] 



"^fe , (142) 



By taking the limit K — > oo, we again see that the low 
frequency behavior of the dynamic structure factor of an 
interacting liquid of hard rods is identical to the oj be- 
havior of the corresponding noninteracting system, pro- 
vided the bare diffusion constant D is replaced by the 
renormalized quantity D, and that the static structure 
of the liquid is taken into account through the factor 
[1 + p.9o(q)]- For finite K, the only effect of the inter- 
nal fluctuations of the vortices on the structure factor on 
long time scales is to introduce the additional "Debye- 
Waller" factor exp(—Tq 2 /y/7CjI). In the limit of non- 
interacting flux lines, where fi — > and D — > Do, the 
exponent (Tq 2 / 'y/K/j,) in this last factor goes to its up- 
per bound (LTq 2 /2AK), the pair distribution function 
go(q) — > S(q), and we recover the result Ijl38|) of an ideal 
gas of vortices. 

Going beyond the above limiting cases, we here would 
like to comment on the experiments of Yao et al. [34| 
and Yoon et al. |35j | who measured the structure factor 
of the vortex liquid in I^S^CaC^Os (BSCCO), and 
attempt a quantitative fit of experimental data using our 
theoretical prediction. For the extremely dilute vortex 
liquid studied in these experiments, the tilt modulus K 
is given by the single vortex value K « eq, which has 
the numerical value 14] e (K/A) = 1.964 x 10 8 /[A(A)] 2 . 
Fig. □ shows a plot of cosh -1 [S(q, z = 0)JS(q, z = L)} 
using the experimental parameters of ref. 35], namely: 



a 
A 
L 
T 



1.5 x 10 4 A, 
0.62 x 10 4 1. 
0.2 mm , 
80K, 



(143a) 
(143b) 
(143c) 
(143d) 



and with the fit parameters a = 1 and r\ = 0.2. Com- 
paring our plot with the experimental curves (Fig. 2 of 
ref. |23), we see that our mean-field approach is able 
to produce a reasonably good qualitative fit of the data, 
which is quite surprising, given the rather simplified form 
of our model Hamiltonian, Eq. IjTjI. and of our analytic 
ansatz for the pair distribution function go(r), Eq. (JISJ. 

It is worth nothing at this point that the confining 
coefficient /i can be reproduced (up to numerical fac- 
tors of order unity) by taking the short wavelength limit 
q — > qsz (short wavelength fluctuations being the dom- 
inant ones in a liquid and qsz = (in/y/3a) being the 
wavevector at the Brillouin zone boundary of a solid at 
the same density) of the compression modulus cn(q) = 
B 2 /47r(l + \ 2 q 2 ) of usual elasticity theory. [3(| That we 
are able to fit the experimental data with a value of the 
compression modulus that qualitatively agrees with elas- 
ticity theory is rather reassuring, and strongly supports 
our claim that an approach based on conformation vari- 
ables of vortices is more adequate to describe flux-line 



liquids than boson mapping methods which by contrast, 
when used to fit the data of Yoon et al. [23 , give a result 
for en which is smaller than the theoretical prediction by 
three orders of magnitude. [IE H3 




FIG. 1: (Color online) 
cosh" 1 [S(q,0)/S(q,L)] vs. q, 
static structure factor, Eq. I13UH . 
experimental curve, from ref. |3a 



Solid line: Plot of 
using our result for the 
Dashed line: approximate 

. As in this last reference, 



the wavevector q is measured in units of qsz = (4tt/\/3o). 

The failure of hydrodynamic approaches to describe 
experimental data in a way that is consistent with stan- 
dard elasticity theory is further indication of the impor- 
tance of separating internal modes and cm. variables of 
the vortices in the liquid state. Indeed, there is a very im- 
portant qualitative distinction between the cm. mode on 
one hand and the internal modes on the other, for while 
the cm. mode is diffusive, internal modes of continu- 
ous media are non-diffusive (due to the elastic restoring 
forces, and this regardless of whether the confining mass 
fi is zero or not). The results of the present study out- 
line the importance of separating the non-diffusive modes 
from the diffusive ones that can be studied using hydro- 
dynamic treatments. 

Another quantity which deserves attention is the 
friction force experienced by the flux liquid driven 
in presence of disorder. For a flux liquid in pres- 
ence of point disorder with correlations A(r, z) — 
Aq exp(— r 2 /2£ i 2 )S(z) /2tt^ 2 we obtain, in the large v limit 
and in d = 3 dimensions: 



7r 



Vt/k i 



(144) 



where factors of order unity have been dropped, and 
where we defined the length £ such that £ 2 = £ 2 + (u 2 ). 
We see that the friction force vanishes at large drives in 
agreement with previous predictions for vortex lattices 
and liquids in the plastic regime near the melting tem- 
perature. 0] 

In summary, in this paper, we have extended the ap- 
proach developed in refs. to the case of a flux liq- 
uid in the presence of a random pinning potential. This 
approach, which is based on the use of the conforma- 
tion variables {r n (z)} as the true dynamical variables in 
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terms of which a Gaussian approximation is taken, gives 
physically more reasonable results 0, |nj than the bo- 
son mapping Q or other hydrodynamic approaches Q 
which, instead, use the density as the basic dynamical 
variable of the vortex liquid. Within our approach, we 
find that the only effect of the pinning potential in the 
static equilibrium case is to renormalize the tilt mod- 
ulus and the confining potential of the internal modes 
of the flux lines, increasing their stiffness and reducing 
their thermal wandering. In a similar fashion, we find 
that in presence of pinning, apart from the appearance 
of nonlinear KPZ terms and standard renormalization of 
the coefficients, the equation of motion of flux lines keeps 
the same form as in an unpinned vortex liquid. As a con- 
sequence, and unlike the hydrodynamic approximations, 
IUlI El we find that the structure factor S(r, z) has 
the same functional form as in a liquid of interacting but 
unpinned flux lines, with suitably renormalized param- 
eters. Our formulation of the equilibrium dynamics of 
vortex liquids is in full agreement with the standard dy- 
namical theory of classical liquids, and through the pair 
distribution function go(r) of the cm. mode, takes into 
account nontrivial correlations in the positions of flux 
lines. In particular, we find that the long time dynamics 
of a liquid of interacting flux lines is qualitatively similar 
to the dynamics of an ordinary classical liquid of hard 
rods with a renormalized diffusion coefficient D (which 
is reduced by the interactions with respect to the free 
value -Do), the only effect of internal fluctuations of flux 
lines at long time scales being to reduce the structure 
factor of the vortex liquid through the introduction of a 
Debye- Waller-like thermal smearing term. 
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APPENDIX A: EFFECTIVE HAMILTONIAN OF 
THE INTERNAL MODES OF FLUX LINES 



shall write: |38| 

<(*)-<oo~(z- z')d z Ui(z) 



(Al) 



upon which one can see that the first term on the right 
hand side of equation l|16f) can be written in the form 



p N 

EE 

a—l i—l 



dz l -8K (d z u?(z)f 



(A2) 



where the long-wavelength disorder contribution SK to 
the tilt modulus of the flux lines is given by 



SK = -- 



i r°° 

— J dzz 2 ViA(r,2) 



r=0 



(A3) 



and vanishes for an ultralocal (in z) disorder, with a vari- 
ance A(r, z) cx S(z). 

Now, for the second and third terms of we need 
to evaluate averages of the form 



(O = </(^-4)> 



(A4) 



where /(r) is an arbitrary function of the space variable 
r. We have 



a,b\ 



Z- 1 Y[dv c m f(r a 0l ~r b 0l )e-^ (A5) 



where Z = Tr( exp(— /3-ff)) . Separating out the integra- 
tions over rj^ and r^ -, we obtain 

</S b }o = / «/(r§i-i)x 

x Z- 1 f H'dr c ok e-^ (A6) 



where the prime in on the product indicates that the 
variables and do not appear in the integration 
measure. Now, two cases have to be distinguished : 

(i) If a = b, i.e. and rL belong to the same replica. 
Then the quantity in the second line on the right hand 
side of equation (|A6|) is given by (see also ref. 



P 



N(N — 1) 



5o (ro, 



L 0j) 



(A7) 



where ga(r) is the pair distribution function of the cm. 
of a given replica in the flux liquid. This leads to the 
following result for (f°f)o 



In this Appendix, we show details of how we perform 
the average (H^l)o, where H^ g is the disorder part of 
the Hamiltonian of the internal modes of the flux lines, 
and where the average is taken with statistical weight 
exp(-SW /T)/Tr(exp-BW /T), As can be seen from 

equation l|16|) . the first term of the Hamiltonian 
does not depend on the cm. coordinates. Assuming 
long wavelength elastic distorsions of the flux lines, we 



N -1 



dr /(r)fifo(r) 



(A8) 



(ii) If a 7^ 6 (rgj and r^- belong to different replicas), 
then the quantity on the second line of equation l|A6(l is 
now given by 



P 

N 2 



4) 



(A9) 
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where go (r) is the pair distribution function of particles 
belonging to different replicas, equation i|27[l . This leads 
to the result 



a b 



>0 = nJ dr f^o(r) 



(A10) 



Using the results (jA8|) - (|A10|) to take the average of the 
second and third term on the right hand side of equation 
(|16J) . and rearraging the resulting sums, we obtain, after 
a few manipulations, the result l|23|) of the text. 



APPENDIX B: HARTREE APPROXIMATION 
FOR THE EQUILIBRIUM DYNAMICS OF 
UNDISORDERED VORTEX LIQUIDS 

In this Appendix, we show how the result of Sec. [V] 
for the effective "mass" term /i, where we used a sim- 
ple Taylor expansion in the flux-linc displacements {u}, 
can be generalized to take into account the potentially 
large excursions of the vortices around their cm. posi- 
tions which are possible in a vortex liquid. Here we shall 
use a Hartree approximation, which is analogous to the 
approach introduced a long time ago in the context of 
spin glasses |3il . and was extended more recently to the 
spherical p-spin model |4jj an d to fluctuating manifolds 
in random media |4lj| . We begin by rewriting the inter- 
action part Ai n t of Eq. (|54|l in the form 



A; 



A 



(o) 



A 



(u) 



where 



AfJ t = Y. J dzdt J V(<l)i<l a iP 0ia (t)e^ 



(Bl) 



(i)-roj(t)] 



x giq-Iuitz.O-Uj'CM)] 



(B2) 



is the cm. part of the interacting action, and where 

i^j q 

(B3) 



Mi" 17 q 

x e »q-K( z >*)-uj(z>*)] 



is the internal modes contribution to Aint- In the cm. 
piece, we get rid of the {u} dependence by making the 
replacement 



Q iq-[Ui(z,t)-u,-(z,t)] 



■ 9q9/3</>° p (z,t;z,t) 



(B4) 



where we 

<t>fAz,t;z',t') 



defined the 

= ([Ui a (z, t) 



correlation function 

Uj a (z',t')][ui/3(z,t) - 



jp{z' , t')]). This yields, for the cm. part A 



(0) 

/III 



•4°t -Y.J dt J V(q)iq a iLf oia (t) e 

i^j q 



iq-[r oi (t)-r j(t)] 



with the effective interaction potential per unit length 
between vortices 



U(q) = V{ci)e-^ q - q ^f (z ^ z ' t) 



(B6) 



which is smeared with respect to the original potential 
V(q) by th thermal fluctuations of the internal modes of 
flux lines. 



We now turn our attention to the relatively more in- 
volved task of constructing a self consistent Gaussian ap- 
proximation for Ant . Following Scheidl (who derived 
a similar self consistent approach for a flux line lattice 
pinned by disorder ^2)' we expand the exponential with 
respect to the displacements and contract the fields in all 
possible ways pairwise until one or two fields remain un- 
contracted. For even and odd terms in the displacement, 
this yields: 



if ia {z,t) x J^y i ic l - MM) - Uj(z,<)]} 2n -> if i c.(z,t) — {~-q a q l3 (j)° :l ll(z,t- 1 z,t)} n 

if ia (z,t) x — -j- ^, {>q- Mz,t) - Uj(M)]} 2 ™ +1 ->■ w>(z,f){iq- [ui(z,t) - Uj(z,t)]} ^{--q a q (t>ij (z,t; z,t)} n 
Resummation yields: 

Ant -E/ dzdtirf(z,t){V«(z 7 t-z,t) + lu^z,t)-u J p(z,t)}V l a f(z,t;z > t)} (B7) 

Mi 
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where we defined for convenience: 

VfJz,t;z,t) = I ^IqJ^e^^W-^W'e-^"'^-^ 2 * 2 '' 1 



V™f(z,t;z,t) = / y(q)^ (?Q %e ^< ^•[ ro •( t )- ro ^ ^ )]e-^' ^ ^' 3 <? (^ • ^;^ ' ^) 

I 



(B8a) 
(B8b) 



The term (t) represents a force acting on particle i 
and arising from particle j, which is proportional to the 
small displacement u". 

With hindsight from the results of Section we shall 
assume that (ui a (z, t)ujp(z, t)) = SijS a /3 (uf a ), which im- 
plies that, for i ^ j , 



1,3 



(z,t;z,t) = —5 a p{u 2 ) 
a± 



(B9) 



The effective smeared interaction potential of Eq. I|B6|) 

becomes 



(BIO) 



The quantities V$ and V°f in Eqs. (B5a ^ -(B55 |l . on the 
other hand, are now given by : 



V&(z,t;z,t) 



V«f(z,t;z,t) 



V(q)iq a e 1 n ' 3 " 2d ± y 1 

2 2 

t/V \- ■ ia-\roi(t)— ro 7 (i)l — — (u ) 

V (q)«9a«9/3e 2d± 



Taking the average of Ai n t over the center of mass posi- 
tions with statistical weight e~ A ° , we obtain the following 
effective interaction action: 

A< ^ ~ T,J dzdt ™?( z > *) {V% {z, t; z, t) + 

+ [u lP (z,t)~u^{z,t)}V^f(z,t;z,t)} (B12) 



where now 

V£i(z,t;z,t) 



(B13a) 

V£/(z,t;z,t) = J V(q)iq a i qf3g o(q)e-^I {u2) 

(B13b) 

where we used the fact that 

lq -[r ,(t)-r 0j (t)]\ _ P_ 



(N-l) 



ffo (q) (B14) 



It is not difficult to see that V a vanishes due to spheri- 
cal symmetry of the interaction potential, V(q) = V(q). 
Putting together all terms in Eq. (|B12|) , one obtains that 
the effective action for the internal modes of flux lines can 



be again rewritten in the quadratic form of Eqs. 1721) - 
(|74b(l . with now the mass coefficient /i given by the self 
consistent equation 



(B15) 



which is identical to the result obtained in ref. [T(i| using 
a static variational approach. 



APPENDIX C: GAUSSIAN HYDRODYNAMICS 
OF FLUX LINES REVISITED 

In this Section, we revisit the Gaussian hydronamic 
formulation of the statistical mechanics of flux liquids, 
which will make it easier for us to compare the results of 
our microscopic approach to the results of previous pub- 
lications 

0111111, which were mostly based on macro- 
scopic coarse-graining methods. We shall first consider 
the case of a liquid of rigid flux lines, before considering 
the general case of a liquid of flexible vortices. 



1. Hydrodynamics of rigid flux lines 

Let us consider a system of rigid flux lines, described 
by the following Hamiltonian: 

H ^X^ofa-r^ + XVe-tfa). (CI) 



where Vo(r) = LV(r), and where we introduced a one- 
body external potential V ex t{r)- The hydrodynamics of 
the liquid described by the above Hamiltonian is con- 
structed in a standard way 0] as follows. First, one 
introduces a variational Hamiltonian, consisting of the 
sum of non-interacting single- vortex Hamiltonians: 



N 



(C2) 



»=i 



with the Hamiltonian Hi to be determined by minimiza- 
tion of the variational free energy F p (the significance 
of the subscript "p" will become clear shortly) , which is 
given by: 



F n 



-ThxZ,, 



(H - H„ 



(C3) 



19 



In the above expression, (•••)« denotes averaging with 
statistical weight e~ l3Hv / Z v , and Z v is the partition func- 
tion associated with the Hamiltonian H v , and is given by: 



Z v = J dri — drva-PX^W, 



dv e-^W 



N 



In the same way, it is easy to show that: 
NfdrHt^e-f 311 ^ 



\H V ) V — 



Jdce-PBi(r) 



(C4) 



(C5) 



Now, the variational average of the density operator 



/?( r ) = (Ei=i <K r - r i))v is given by: 

Ne -PHi(r) 



p{v) 



so that 



in (4?) = ^ i(r)_in (/ dr 



-0fli(r) 



(C6) 



(C7) 



Hence: 



- NT\n( f dre- pHl{r) 



(C8) 

Comparing the above equation with Eqs. I|C4|I and l|C5() . 

we see that: 

j dv Tp(r) In (^) = -TlnZ„ -<#„)„. (C9) 

On the other hand, it is not difficult to show that the 
variational average of the original Hamiltonian H is given 
by: 

(H) v = 1 /" drdr'p(r)F (r-r')p(r') 

+ /" dr Ve»t(r)p(r) . (CIO) 



Collecting together the results (|C9(1 and (|C10(1 , it follows 
that the variational free energy of Eq. (|C3(I can be writ- 
ten as a functional of the averaged density p(v) (hence 
the subscript "p" in F p ), and is given by: 

F P = J drTp(r)ln(^) + J dv V ext {v) p{v) 

+ ^[ dvdv' p(v)V (v-v')p(v'). (Cll) 

The first term on the right hand side of the above equa- 
tion is the entropic contribution to the free energy of the 



liquid, which is qualitatively important (we shall see be- 
low that this term modifies the temperature dependence 
of the density response and correlation functions — see 
Eqs. (|C18|1 and (|C20(l ). but which has been systemat- 
ically overlooked in previous studies of flux-line liquids. 
Li LJ 13 We now need to minimize the variational free 
energy F p with respect to the density p(v), which may 
be thought of as an "order parameter" , with the con- 
straint that the total number of particles N is held fixed. 
This amounts to minimizing the variational version of the 
grand potential: 

Q p =F p -p J dvp(v), (C12) 
with respect to variations in the density p(v), using [i^| : 



0, 



(C13) 



with the Lagrange multiplier p (chemical potential) in 
Eq. ljUT2)l fixed by the condition / dv p(v) = N. This 
minimization procedure leads |43j to an expression for 
the one-body density p(v) which is identical to Eq. (|C6|) 
above, with H\ given by: 

H^v) = V ext (v) + J dv' V Q (v - r')p(r') . (C14) 



Now, the density response function is given by: [4j 

Sp{v) 



X(r,r') 



6V ext (v>) 



(C15) 



Taking the functional derivative of p(v) in Eq. (|C6|) with 
respect to V ex t(v'), one can easily show that: 



X(r,r') = (3p(v) 



5H x {v) 



5V ext (v>) 

(3 P (v) r . ,5zri(n) 



y dv! p(ri)- 



N J Lry J 5V ext (v') 
SH^v) 



pp(v) 



6V ext (v>) 



(C16) 



where, in going from the first to the second equality, we 
discarded a term proportional to 1 /N which vanishes in 
the thermodynamic N — > oo limit. Now, using the result 
(fC14l for Hi into Eq. fC16jl . we obtain: 



X (v, r') = (3p(v) [5(v - r') - / dv" V (v - r")x(r", r')] . 

(C17) 

In the absence of an external one-body potential (i.e. 
when V ex t = 0), the system is translationally invariant, 
and x(r, r') = x( r — r ')- Taking the Fourier transform of 
Eq. (|C17(I then leads to the result: 



X(q) 



1 



V (q)+T/p 



(C18) 
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The T/p term in the denominator originates from the 
p\a.p entropic term in expression (|(Jllfl of the free energy. 
As we mentioned above, this term, which determines the 
temperature variation of the density response function, 
was totally ignored in previous studies of flux liquids. 
Pi El til From the above response function, the Ursell 



function S nn {v) = (p(r 
and notation of ref. 41 



)(0)) —p 2 (we use the terminology 



|) is given by: 
S nn {v) = T X (r) , 
which gives, in Fourier space: 

S n „(q)=T X (q) 



V (q)+ T/p 



(C19) 



(C20) 



The above Ursell function can be obtained directly from 
Eq. (|C11|I by writing the density p(r) as p(r) = p + 
Sp(r), and expanding F p up to quadratic order in the 
density fluctuation Sp(r). Taking the Gaussian average 
of the product 5p(r)5p(0) by integrating directly over the 
density fluctuation field Sp(r), 



(6p(r,z)5p(O,0)) 



f[d(5p(r, z))}Sp(r,z)Sp(O,0) 



-0F P 



f[d(6p(v,z))}e-^ 

(C21) 

leads directly to the result IJC20I) . Note however that 
this last procedure to obtain the Ursell function is an 
ad hoc one (that is valid only because f2 p is stationary 
with respect to variations of the density p(r) around its 
equilibrium value p), and that the most systematic and 
justifiable way to obtain S nn in this variational formula- 
tion of hydrodynamics is through extremizing the grand 
potential fl p with respect to the density (a step that 
is rigourosouly exact |4J|), and then using the resulting 
variational Hamiltonian to find the density response func- 
tion and hence S nn (q), as we did in Eqs. HC11|) - HC20(1 . 
We insist that in the latter method no integration over 
density variables (as in Eq. (|C2f() ') is performed, and we 
thus see that Eq. (|C2f() is by no means a necessary step 
to obtain the Ursell function. 



2. Hydrodynamics of flexible flux lines 

We now want to generalize the approach of the above 
Subsection to flexible flux lines. For the sake of homo- 
geneity with the rest of the paper, we shall use the same 
Hamiltonian as in the text, namely: 

N r i <7 2 

H = E I dz k K \~dj) + v ^(^)^)] 

+ iE/ dzViTiW-Tjiz)), (C22) 

where we again, by analogy with the previous Subsection, 
have introduced a one-body external potential V ext (r, z). 
In an obvious generalization of Eq. (|C2|) to continuous 



systems, we shall use the following variational Hamilto- 
nian: 

N 

H v = J2Hi[ri(z)}, (C23) 

i=l 

where Hi now is a functional of the entire trajectory rj(z) 
of the i-th flux line. The corresponding partition function 
is given by: 

r N 

Z v = J [dr ! (z)] • ■ ■ [dr N (z)} exp(-/3EtfiM*)] 



[dr(z)] e-^I'Wl 



/V 



(C24) 



The variational average of H v on the other hand is given 
by: 

_ Nj[dr(z)}H 1 [r(z)}e-^^ 
Wv)v- f[dr(z)] e-MM*)) ' 1 ' 

in total analogy with Eq. (|C5p . Let us now define the 
operator: 

N 

p[r( Z )]=^n%W^,W), (C26) 

i=l z 

whose statistical average (p[r(z))} may be interpreted as 
the average probability density for finding a vortex with 
a specific trajectory r(z) in the system. Now, the vari- 
ational average p[r(z)] = (p[r(z)]) v of this trajectory- 
density operator is given by: 

N -PHlM*)] 

p ^ z)] - jwm e-^*a ■ (C27) 

Forming Tp[r(z)]\n(p[r(z)]/N), and then taking the 
functional integral over the trajectory r(z), one can again 
easily show that: 

[dr(z)} Tp[r(z)} In (^^) = -TlnZ v - (H v ) v . 

(C28) 

Now, the variational average of H can be expressed in 
terms of the trajectory density p[r(z)], and is given by: 

+ V ext (r(z),z)]p[r(z)} 

+ p[r(z)}V(r{z)-r'(z))p[r'(z)]}, (C29) 

and hence we obtain that the variational free energy F p — 
—T\iiZ v + (H — H v ) v can be written in the form: 

'p[r(z)}' 



(H) 



F p = J dz{ J [dr{z)\ Tp[v(z)]\n 



N 





'dr\ 








-.dzJ 



[dr(z)}[dr'(z)]p[v(z)}V(v(z) - r'(z))p[r'(z)] }, 

(C30) 
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and is a generalized functional of the averaged trajec- 
tory density p[r(z)]. In the subsequent steps of the hy- 
drodynamic method, which involve minimizing F p with 
respect to the density p[r(z)], it is convenient to dis- 
cretize the z axis, with unit step d, and think of p[r(z)] 
as a function p(r(Q),r(d), . . . ,r(Md)) of the positions 
(r(0), r(d), . . . , r(Md)) of M+l particles interacting with 
each other harmonically, with a spring constant K/d. 
The integration measure in this case can be defined as: 



M 



[dr(z)] = Y[ dr{md) . 
Now, from Eq. I|C27|I we see that: 



[dr(z)} p[v{z)]=N, 



(C31) 



(C32) 



and hence we see that the grand potential il p = F p — pN 
can be written in the form: 



F p -p [dr(z)] p[r(z)] 



(C33) 



In the same way as in the previous Subsection, 
extremization of Q p with respect to p[r(z)] — 



p(r(0), r(d), . . . , r(Md)) (thought of as a "generalized" 
order parameter) under the constraint l|C32|) leads to the 
following result for the variational Hamiltonian Hi[r(z)]: 

HMz)} = Jdz{[lK(^y + V ex Mz)) 

+ J [dr'(z)} V(r(z) r'(z))p[r'(z)]}.(C34) 

Like we did in the text, we now decompose the trajectory 
r(z) into cm. and internal coordinates, r(z) = ro + u(z), 

with ro = J Q dz r(z) being the cm. coordinate of the 
trajectory r(z), and Taylor expand the interaction part 
in the small displacement u(r): 

V(v{z)~v'{z)) = V(r -r'(z)) 

+ u a (z)d a V(r - r'(z)) 

+ ^-u a (z)up(z)d a d l3 V(r () - r'(zj) . 



(C35) 



The Hamiltonian Hi becomes: 



//,[..'(:)] = / dz[~K(^) 2 + V ext (r(z))] + J dz j[dv'(z)\ V(r - r> (z)) p[r' (z)} + 

dzu a {z) [[dv'(z)}d a V(r -r'(z))p[r'(z)] + ~ J dz u a (z)u (z) J [dv'{z)]d a d p V{y Q -v'{z))p[v'{z)). 

(C36) 



Using the fact that J dzu(z) = 0, we see that the third As defined above, the quantities p, a p depend on ro. We 
term on the right hand side of the above equation van- shall however verify a posteriori that in the homogeneous 
ishes, and therefore that Hi[r(z)] in the absence of an liquid state this dependence drops out and the Pap's re- 
external potential (V ex t = 0) can be written in the form: duce to ordinary constants (and are in fact all equal to 

zero in the hydrodynamic limit). 

An immediate consequence of the decomposition l|C37|l 
is that the density matrix also can be written in a decou- 
w h ere . pled form: 



Hi[t(z)] = ffi,c.m.(r ) + H hu [u(z)] , 



(C37) 



ffi,cm.(r ) = / dz I [dr'(z)]V(r - r'(z)) p[r'(z)} , 

(C38a) with: 



p[r{z)] = Pernio) p u [u(z)] 



H hu [n( Z )] = J dz 



-p a pu a (z)up(z) 



(C38b) 



In the above expression of Hi u [u(z)] , we defined the con- 
fining "mass" tensor: 



Pc.m.( r o) — 
Pu[u(z)} = 



-/3ffl,c. m .(r ) 



Jdr e -/SHi,..».(ro) 

e -^l,„[ U (2r)] 



(C40) 

(C41) 
(C42) 



J[du(z)} e -PHiM»)] ' 

fu 't ma q t / 1 \\ r w , PQn , Now, if we repl ace th e density matrix p[r(z )] by the 

p aP = J [dr (z)} d a d V{r o - r (z))p[r (z)} . (C39) decoupled form ^ back into expression JH381 of 
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Hi tC , m ,, we obtain: 

ffi, c .m.(ro) = Jdr' LV{r -T' ) Pc . m Xr' ), (C43) 
where we defined: 

U(r) = /"[du'(z)] y(r-u'( 2 )) /9u [u / ( 2 )],(C44a) 
= y e ik - r -^ 2<u2 V(k), (C44b) 

and where, in going from the first to the second line of 
the last equation, we assumed that the density matrix 
p u [u'(z)] represents an isotropic Gaussian distribution 
for the displacement field {u'(z)}. In Eq. l[C44aJI . the 
integration measure [du(z)] stands for: 

M-l 

[du(z)] = Y[ du(md) . (C45) 

m— 

Note that the transformation from the variables {r(md)} 
(to = 0, ...,M) to the variables {ro,u(md)} (to = 
0, . . . , M — 1) being linear, it has a constant Jacobian, 
that we shall henceforth ignore for simplicity. 

Now, if we insert the decoupled form (|C40(I of /9[r(z)] 
into expression (|G39|) . we obtain: 

Pa/3 = J dr' Q p c .,nX r o) 

x J[du'(z)\ d a df)V(r -r' - u' (z)) p u [u' (z)] , 

= f dr' o d a d V(r o -r' o )p c . m Xr' o ). (C46) 

In a uniform flux liquid, the averaged cm. density re- 
duces to a constant, p c .mX r 'o) = Pi an d the above in- 
tegral vanishes. This is an direct consequence of the 
coarse-graining procedure, for if instead of the averaged 
cm. density p c .mX r 'o)-> we were still dealing with the cm. 
densit y operator p c . m X r o) = J2iLi ^( r o - r 0i), then Eq. 
(|G46I) would give: 

N 

AV3 = ^nr -r 0j ), (C47) 

i—l 

which is very similar to the undisordered version of Eq. 
(|12fl . which then would yield finite and isotropic "mass" 
coefficients [x a p = p5 a 0. 

Now, if we calculate the structure factor S(r, z) by di- 
rectly taking the average of (p(r, z)p(0, 0)) as in Eq. l(F7|) 
and using the decoupled density matrix of Eq. (|C40(I , one 
can easily show that the structure factor 5(q, z) of the 
flux-line liquid has the form given in Eqs. H127fl and 113U|) 
of the text, with /i = 0, namely: 

S(q, z) = pigoid) e-f^' d - + pexp ( - ^\z\) . 

(C48) 



We have: 

p 2 5o (q) = (2^)V <5(q) + S ,„„(q) - p , (C49) 

where So lV , n is the Ursell function of the cm. mode 
(S , nn {r) = (pc. m {r)Pc. m X )) - p 2 ) , which can b e eas - 
ily derived from the cm. effective Hamiltonian (|C43(1 . 
following the same steps as in the previous Subsection, 
with the result: 

s °" (q) = Z^7/ (C50) 

Using the result (|C49|) . Eq. (|C48|I can be rewritten in 
the form: 

S(q, z) = (27r) 2 p 2 S(q) + [5 ,„„(q) - p] e^^' d ^ 

+ pexp(-|^|z|). (C51) 

Since p, = in the hydrodynamic limit, the mean squared 
displacement (u 2 ) oc L, and hence e~ q ^ d± is expo- 
nentially small for practically all values of q such that 
< q < 1/a. Hence the above expression of the struc- 
ture factor becomes: 

S(q, z) ~ (2t:) 2 p 2 5(q) + pexp ( - ^\z\) , (C52) 

which implies that the Ursell function of the vortex liquid 
is given by: 

5'nn(q I 2) = pexp^-^|z|^, (C53) 

and is identical, in the hydrodynamic limit, with the 
Ursell function of an ideal gas of non-interacting flux 
lines. We thus see that, while the present formulation of 
Gaussian hydrodynamics yields a structure factor that is 
of the correct general functional form, Eq. (|G48(I (unlike 
previous formulations of refs. 0,0,0,0); it fails to pro- 
duce a non-zero value for the confining "mass" term, with 
the consequence that the ensuing strong fluctuations of 
flux lines completely smear out the effect of interactions 
between flux lines (the second term on the right hand 
side of Eq. ESQ )- 

The decomposition given in Eq. I|C23|) . which approx- 
imates the TV-body Hamiltonian of the system by a sum 
of N noninteracting one-body Hamiltonians, is not the 
only possible choice for the variational Hamiltonian H v of 
the hydrodynamic method. One other (but not necessar- 
ily equivalent) possible choice consists, in the discretized 
scheme where the z axis is cut into M + 1 equidistant 
slices, in using the ansatz: 

N M 

h v = h ( r ( md )> md ) . ( C54 ) 

i—l m— 

with the Hamiltonian of a single vortex "element" 
h(r(md),md) at height z = md to be determined varia- 
tionally. The above ansatz amounts to assuming for the 
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Hamiltonian 7Ji[r(z)] of the preceding paragraphs the 
following form: 



M 



iJi[r(z)] = h(r(md),md) 



(C55) 



m=0 

The variational partition function Z v is now given by: 

M 



Z v =f[(fdr, 



- Ph(r,rnd) 



N 



(C56) 



while the variational average (H v ) v is given by: 



x ' jdr e -l3h(r,md) ^ ' 



If we define the density operator at height z — md, 
p(r,md), by: 



JV 



){r,md) = ^<5(r-r i (md)), (C58) 



then it follows that the averaged density p(r, md) 
(p(r, md)) v is given by: 



p(r, md) = 



jy e -0h{r,md) 



J dr e -PH r > rnd ) 
Thus, here again we can write: 



(C59) 



dr Tp(r, md) In ^- 



p(r, md) \ N J dr h(r, md)e- f3h ^ T ' md '> 



N 



J dr e -Ph(r,md) 
I 



TN In / dr e 



-ph(r.md) 



(C60) 



Summing over m, we obtain 

M 



Y dr Tp(r, md) In ( 

m=0 



p(r, md) 
N 



(C61) 

Taking the variational average (-ff% of the original 
Hamiltonian H, we finally find that the variational free 
energy F p = —Tin Z v + (H — H v ) v can be written in the 
form (we now switch back to a continuum notation): 



F p = I drdz 



T 



d p { r,z)\n(^l) + 1 -K 1 t\r,z) 

+ 1 ~Jdrdr>fdz P (r,z)V(r-r>) P (r',z) 

+ J drdzV ext (r,z)p(r,z), (C62) 

where K\ = K/p, and where we defined the "tilt" field 
operator 0,3: 



JY 



dr.: 



t(r,*) = £^«(r -*(*)) 



(C63) 



The tilting field t(r,z) and the density field p(r,z) are 
not independent, but_ are_related to each other by the 
continuity equation 



d z p(r, z) + V_l • t(r, z) = . 



(C64) 



Nelson and coworkers 0, Q obtain density and tilt cor- 



relation functions by expanding F p of Eq. I|C62|) to 
quadratic order in Sp(r, z) = p(r, z) — p (omitting, for 
some unstated reason, the entropic p In p term) , and cal- 
culating statistical averages in the manner ofEq. ltC2TTl . 
with the constraint 1C64|) enforced. For example, for the 
Ursell function S nn (r,z) — (Sp(r, z)Sp(0, 0)), these au- 
thors write: 



S nn (r,z)) = 



J[d(6p(r, z))] /[rft(r, z)] 5(d z p(r, z) + V ± - t(r, z)) 5p(r, z)8p{0, 0) 
J[d(6p(r, z))} J[dt(r, z)] 5(d z p(r, z) + V±- t(r, z))e'^ 

I 



(C65) 



with the result: 



Tql 



(As we mentioned already earlier, in previous treatments 
the T I pd term in the denominator of the above expres- 



[V(d)+T/(pd)]ql+K iq r 



(C66) 
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sion, which comes from the entropic plnp term in F p , is 
missing.) Unfortunately, the Gaussian integration in Eq. 
(IC65II is not justifiable from the point of view of a vari- 
ational approach. In other words, correlation functions 
obtained by using Eq. (|C65|I cannot be reproduced by a 
standard variational method, which here would consist in 
extrcmizing the free energy F p with respect to the fields 
5p(r,z) and t(r, z), with steps similar to those of Eqs. 
(|ClHl - <|C20f) of the previous Subsection. Indeed, as we 
have seen in the end of Subsection IC II above, the ratio- 
nale behind the Gaussian averaging in Eq. I)G65|1 is that 
the free energy is (presumably) a functional of 6p(r, z) 
and t(r, z) (with the constraint (|G64|) enforced) that is 
stationary at thermal equilibrium. 

Technically, the standard way to implement the con- 
straint i|G64|) in a variational procedure is through the in- 
troduction of a (functional) Lagrange multiplier A(r, z), 
whereby one defines the following, modified grand poten- 
tial (the last term in this equation is simply —p,N): 

flp = F p + J drdz \{r,z)(d z p(r,z) + V_ L • t(r, zfj 

- | J drdz p(r, z) . (C67) 

Extremizing fl p with respect to p(r,z) and t(r, z) leads 
to the following coupled equations: 

d z \(r,z) = ^(l n (^^)+l)+U e:ct (r,z) 

+ J dv' U(r-r>(r',z) + |, (C68a) 
V±A(r,z) = #it(r,js). (C68b) 

Eq. (|C68b|l is a vector equation of standard form, which 
amounts to finding the "electric potential" A(r, z) asso- 
ciated with the planar "electric field" t(r, z), and has a 
well-defined solution in A if and only if V j_ x t = 0. Since 
t(r, z) is a randomly fluctuating field that does not nec- 
essarily satisfy this last condition, we arrive at the very 
important conclusion that it is not legitimate to extrem- 
ize the grand potential with respect to the pair of vector 
fields {p(r, z), t(r, z)}, and hence that it is not legitimate 
to calculate statistical averages using the procedure ex- 
amplified in Eq. (|CJ65|) . 

We thus see that the previous formulations of Gaussian 
hydrodynamics which lead to expressions for the Ursell 
function of the form given in Eq. (|C66|1 correspond to 
a convoluted and ad hoc attempt, without any ratio- 
nale other than hand-waving symmetry considerations, 
to generalize the variational Gaussian hydrodynamics of 
point particles to continuous systems. In fact, even if 
we ignore this lack of rationale and accept the use of Eq. 
(|G65ll , one other source of inaccuracy of the previous for- 
mulations of Gaussian hydrodynamics lies the underlying 



Hamiltonians of Eqs. ljG54|) - (|C55(l . which these theories 
are all implicitly based on. These Hamiltonians indeed 
represent a very crude approximation to the Hamiltonian 
of a single flux line in the first place, since Hi[r(z)] in- 
volves relatively strong (harmonic) interactions between 
flux line segments while Eq. I|G55[) models a single flux 
line as a superposition of non-interacting elements. In 
fact, a necessary condition for the applicability of the 
Gaussian hydrodynamic approach is that the interactions 
be weak. |43j For continuous systems with constituent 
parts interacting strongly, a correct formulation of Gaus- 
sian hydrodynamics must take these strong elastic inter- 
actions into account as exactly as possible, as we did in 
Eqs. (|C23|I - <|G34() . for otherwise one may obtain abnor- 
mal behaviour for z correlations in the system (for ex- 
ample, the Ursell function S nn (q,z) obtained from Eq. 
(|G66ll decays more rapidly than the corresponding quan- 
tity for an ideal gas of non-interacting flux lines, which 
is very surprizing, as discussed in detail in ref. flO| ~) . 



3. Gaussian hydrodynamics of flexible flux-lines in 
presence of disorder 

We now generalize the formulation of Gaussian hydro- 
dynamics that we developed in the previous Subsection 
to flexible flux-lines in presence of disorder. Since the 
disordered case involves only minor technical modifica- 
tions of the undisordered hydrodynamics, we shall only 
give the salient features of the calculation, leaving out 
the (obvious) technical details. Our starting point is the 
replicated Hamiltonian of Eq. (jSJ) , which we rewrite here 
for clarity (we remind the reader that p denotes the total 
number of replicas) : 

S = EE HMir) + E vtfw-w)} 

o=l »=1 J j&i) 
P N 

" 2T £ £ Jdz^'HrU^-r^z-z'). 

a, 6—1 i,j=l 

(C69) 

Rewriting the above Hamiltonian in terms of the trajec- 
tory density operator of Eq. JU26J), and using a varia- 
tional ansatz for the total Hamiltonian of the system of 
the form: 

p N 

ff« = ££#i [<(*)], (cro) 

a—l i—1 

it is not difficult to show that the variational free energy 
F p = -T In Z v + (H - H v ) v is given by: 
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F p = ]T / 'dz |[d r ( z )]{ T p>(z)]ln(^^ 

a—1 

+ dz j[dv{z)][dv'(z)]p a [v{z)]V{v{z)-v'{z))p a [v'{z)] 

a—1 

- J dzdz' J[dv{z)][dv\z')]p a [v{z)]^{v{z)-v\z'),z-z')p b [v'{z')]. 



(C71) 



Extremizing F p with respect to p a [r(z)] = p a (r(0), . . . ,r(Md)) leads to an expression similar to the right hand side 
of Eq. lfC27| for p a [r(z)], with the following expression for the effective Hamiltonian Hi: 



ffi [*(*)] = j dz {[l K (^)' * + Ve*Mz))}+ j[dr'(z)]V(r(z)-r'(z))p a [r'(z)] 
- J dz' J[dr'(z')} A(r(z)-r'(z'),z- z')p b [r'(z)}} . 



(C72) 



We now write r(z) = r + u(z) and Taylor expand the above Hamiltonian up to quadratic order in the small 
displacement field u(r, z). Here again we find that, to order 0(u 2 ), Hi can be written in the decoupled form (|C37|) . 
with (here A(r) = dz A(r, z)): 



Hi,cm.(ro) = [dz [[dv'(z)]V(r -r'(z))p a [r'(z)]-^Yl j dz 1 j [dr' {z')]A{v Q - v' {z'),z - z')p h [v' [z')} , 
Hi. u [u(z)] = [ dz \ K {^) + ^Vapu a (z)up(z) 



(C73a) 
(C73b) 



In the above expression of ifi iU [u(z)], the confining "mass" tensor is given by: 

^ = J[dr'(z)}d a d p V(v -v'(z))p a [r'(z)}-^J2 J dz' J [dv'(z')}d a d p A(r - r'(z'), z - z')p b [v'(z')} 



(C74) 



r 



and vanishes in hydrodynamics for the same reason as 
in the undisordered case. Also, like in the pure case, 
the density matrix decouples, p[r(^)] = p c .m. [i"o] jO u [u(z)], 
and hence the internal modes {u'} in the expression of 
Hi^ c .mX T o) can be integrated out, with the result: 



Hi, c . m .(r ) = V e a (r ) + J2 [ dr 'o T ^o - r' )p b . m .(r') , 

b J 

where we introduced an external "source" potential V" 
that depends only on the cm. position ro, and where the 
kernel T a b is given by: 



In the above equation, the potential V is given by Eq. 
(IC44b|l . while A is similarly given by: 

A(r) = [[du'(z)}A(v-u'{z))p a u [u'(z)}, (C77) 



3 lk - r -^ 2 < u2 >A(k). (C78) 



By analogy with Eqs. I|C15(I - (|C17(I . the density response 
function x a6 (r,r') = -5p^ m .(r )/6V b (r' ) satisfies the 
following equation: 

X ab (r ,r' ) = 0p a (r)[5 ab S(r-r') 

- I <r ac (r -r' ') X cb K,r' ) 



: a6 (r ) = L(y(r )5 ab 



A(r ) x 
T J 



(C76) 



(C79) 

For a homogeneous liquid (V e — 0, p(r) = p = Cst.), 
X ab (ro,rQ) is translationally invariant, and the above 
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equation can be cast, in Fourier space, into the following 
matricial form: 



E fac (q)x cb (q) = <u 



with 



~ , s fT \ LA(q) 
r ac (q) = (- + LV(q)) S ac ^ 



(C80) 



(C81) 



The matrix T can easily be inverted using an identity for 
inverting p x p matrices of the form: 



namely: 



(A = aSij + b, 



A - -5 - 

y a y a(a+p6) 



(C82) 



(C83) 



For the response function x a& (cOi this gives in the limit 
p — > the following result: 



X Qb (q) 



1 



LA(q) 



LV(q)+T/p T[iV-(q)+T/p] 



r(C84) 



and hence the diagonal (in replica space) Ursell function 
for the cm. mode is given by: 



So,nn{q) 



T 



LA(q) 



LV(cO+T/p [LV(q)+T/p]' 



(C85) 



We thus see that disorder produces a Lorentzian-squared 
correction to the Ursell function fo the cm. mode in 
the hydrodynamic limit. The Ursell function of the flux- 
line liquid is however unchanged with respect to the pure 
case (since the confining mass p is still zero), and is given 
by Eq. HC53I) . in contrast to the results of refs. [3j and 
Q. Note, however, that if a finite mass coefficient p, then 
the Ursell function of the flux liquid in the hydrodynamic 
limit is given by Eq. I|C51|) , which in the present context 
becomes: 



Snn (qj z ) 



T 



LA(q) 



LV(d) + T/p 

,-9 2 (« 2 >/d 



[LV(cO+T/pY 
Tq 2 



+ pexp(--£\z\). (C86) 
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